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o : 

The origin of the phenomenological deterministic laws that approximately govern the 
quasiclassical domain of familiar experience is considered in the context of the quantum 
mechanics of closed systems such as the universe as a whole. A formulation of quantum 
mechanics is used that predicts probabilities for the individual members of a set of alterna- 
tive coarse-grained histories that decohere, which means that there is negligible quantum 
interference between the individual histories in the set. We investigate the requirements 
for coarse grainings to yield decoherent sets of histories that are quasiclassical, i.e. such 
that the individual histories obey, with high probability, effective classical equations of 
motion interrupted continually by small fluctuations and occasionally by large ones. We 
discuss these requirements generally but study them specifically for coarse grainings of the 
type that follows a distinguished subset of a complete set of variables while ignoring the 
rest. More coarse graining is needed to achieve decoherence than would be suggested by 
naive arguments based on the uncertainty principle. Even coarser graining is required in 
the distinguished variables for them to have the necessary inertia to approach classical 
predictability in the presence of the noise consisting of the fluctuations that typical mech- 
anisms of decoherence produce. We describe the derivation of phenomenological equations 
. £h ! of motion explicitly for a particular class of models. Those models assume configura- 

^ tion space and a fundamental Lagrangian that is the difference between a kinetic energy 

quadratic in the velocities and a potential energy. The distinguished variables are taken 
to be a fixed subset of coordinates of configuration space. The initial density matrix of the 
closed system is assumed to factor into a product of a density matrix in the distinguished 
subset and another in the rest of the coordinates. With these restrictions, we improve the 
derivation from quantum mechanics of the phenomenological equations of motion govern- 
ing a quasiclassical domain in the following respects: Probabilities of the correlations in 
time that define equations of motion are explicitly considered. Fully non-linear cases are 
studied. Methods are exhibited for finding the form of the phenomenological equations 
of motion even when these are only distantly related to those of the fundamental action. 
The demonstration of the connection between quantum-mechanical causality and causality 
in classical phenomenological equations of motion is generalized. The connections among 
decoherence, noise, dissipation, and the amount of coarse graining necessary to achieve 
classical predictability are investigated quantitatively. Routes to removing the restrictions 
on the models in order to deal with more realistic coarse grainings are described. 
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I. Introduction 

In a universe governed at a fundamental level by quantum-mechanical laws, character- 
ized by indeterminacy and distributed probabilities, what is the origin of the phenomeno- 
logical, deterministic laws that approximately govern the quasiclassical domain of everyday 
experience? What features and limitations of these classical laws can be traced to their 
underlying quantum-mechanical origin? This paper addresses such questions in the con- 
text of the quantum mechanics of closed systems - - most realistically and generally the 
universe as a whole. 

It is a familiar notion in the quantum mechanics of simple "measured" systems that 
some coarseness in their description is needed if they are to approximate classical behavior. 
The Heisenberg uncertainty principle, for example, limits the accuracy with which posi- 
tion and momentum can be specified simultaneously. Successions of such suitably crude 
measurements of position and momentum can be correlated by the classical equations of 
motion following from the fundamental action of the system, provided there is a suitable 
initial state — typically a narrow wave packet. 

In discussing the quasiclassical domain of familiar experience, however, we are dealing 
with a much more general situation than is envisioned by elementary analyses of the above 
type. We are concerned, first of all, with the classical features that may be exhibited 
by the behavior of the closed system irrespective of whether those features are receiving 
attention from "observers". We should be able to deal, for example, with the classical 
behavior of the moon whether or not any "observer" is looking at it. Second, we are 
necessarily concerned, not just with classical behavior exhibited by correlations among 
events at a few times of our choosing, but also with the classical behavior of whole orbits 
in as refined a description of the system as is possible. Third, we are concerned with 
phenomenological equations of motion, the form of which may be only very indirectly 
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related to that of the fundamental action. The fundamental action, after all, may be 
that of heterotic superstring theory, while the equations of motion (such as the Navier- 
Stokes equation) governing the familiar quasiclassical domain involve such quantities as 
the averages of the densities of field energy and momenta over volumes very much larger 
than the Planck scale! Further, the applicability of effective classical equations of motion 
may be branch-dependent, that is, contingent on events that have happened. The classical 
equations governing the motion of the moon, for example, are contingent on its actually 
having formed in the early history of the solar system. In such general situations, simply 
identifying the form of the phenomenological classical equations of motion becomes an 
important problem. 

It is a characteristic feature of the general situations described above that a much 
coarser graining is needed for quasiclassical behavior would be naively suggested by ar- 
guments based on the uncertainty principle. As we shall argue below, a large amount of 
coarse graining is needed to accomplish decoherence, which is an important ingredient of 
quasiclassical behavior as well as a sufficient (stronger than necessary) condition for as- 
signing probabilities to the coarse-grained histories of the closed system. Further coarse 
graining is then necessary to achieve the "inertia" required for approximate predictability 
in the presence of the noise from the fluctuations that typical mechanisms of decoherence 
involve. All this coarse graining has important consequences for the form of the effective 
classical equations of motion. Their form may be as much influenced by the character of 
the coarse graining and the mechanisms of decoherence as by the fundamental equations of 
motion. The effective classical equations of motion necessarily include phenomena like dis- 
sipation arising from the mechanisms that produce decoherence. This paper is concerned 
with general methods of deriving the form of the phenomenological classical equations of 
motion and with the description of the noise that causes deviations from those equations 
and from classical predictability. As a result of the coarse graining, the noise includes the 
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effects of classical (typically statistical-mechanical) fluctuations as well as quantum fluctu- 
ations, and these effects are mixed. The resulting indeterminacy, as indicated above, goes 
far beyond the elementary indeterminacy of the Heisenberg uncertainty principle. An ac- 
curate framework for prediction may be achieved by incorporating a generalized Langevin 
force that represents this noise into the classical equation of motion. Our paper is therefore 
concerned with the derivation of the general form and distribution of such forces as well 
as with the equation itself. 

It is known, of course, that even in the classical deterministic limit one can encounter in 
non-linear systems the phenomenon of chaos, in which the sensitivity of the outcome to the 
initial conditions is exponential in time. In the presence of chaos, even small fluctuations 
(including quantum fluctuations) can be amplified to produce large uncertainties in later 
behavior. A treatment of this combined effect of classical chaos and of indeterminacy 
arising from quantum mechanics, including the associated coarse graining, has often proved 
elusive in discussions of quantum chaos, but is amenable to analysis by the methods we 
shall describe [1]. 

II. Decoherence, Inertia, and Equations of Motion 

In this section we give a qualitative discussion of the role of decoherence in the deriva- 
tion of phenomenological classical equations of motion. This will serve to review some 
aspects of the quantum mechanics of closed systems, motivate the subsequent mathemat- 
ical derivations, and make connections with earlier work known to us. For the quantum 
mechanics of closed systems, we follow our discussion in [2,33,3], where references to the 
earlier literature may be found. 

Most generally, quantum mechanics predicts the probabilities of the individual mem- 
bers of a set of alternative, coarse-grained, time histories of a closed system. By a coarse- 
grained history we mean, for example, one for which not every variable is specified and 

4 



those that are specified are not fixed at every time or with arbitrary precision. It is evident 
why coarse-grained histories are of interest to us as observers of the universe. Our obser- 
vations fix only a tiny fraction of the variables describing the universe and fix those only 
very imprecisely. As observers we therefore necessarily deal with a very coarse-grained 
description of the closed system in which we live. However, from the theoretical point 
of view, it is not necessary that the description be so very coarse-grained or that the 
coarse graining be so dependent on us. There is a more fundamental reason for interest in 
sets of coarse-grained alternative histories: In the quantum mechanics of closed systems, 
probabilities may be assigned only to those sets of histories for which there is negligible 
quantum mechanical interference between the individual histories in the set (given the 
system's Hamiltonian and initial quantum state) [5, 6, 2]. We shall define various kinds of 
decoherence [33], all of which imply the vanishing of this interference and for all of which 
coarse graining is necessary. An arbitrarily fine-grained description of the universe would 
reveal the phase correlations between histories, while in a coarse-grained description they 
may be absent. The probabilities of decoherent sets of coarse-grained alternative histories 



The term "decoherence" is used in several different ways in the literature. We have used the 
term to refer to a property of a set of alternative histories of a closed system. Specifically 
a set of coarse-grained histories decoheres when there is negligible interference between 
the individual histories in the set as measured by one of the several conditions discussed 
in this Section. In the literature the term "decoherence" has also been used to refer to 
the decay in time of the off-diagonal elements of a reduced density matrix defined by a 
coarse-graining of variables at a single moment of time, for example, the density matrix 
defined by (4.21). These two notions are not exactly the same. A reduced density matrix 
can be defined for those special coarse grainings that distinguish a fixed set of coordinates. 
However, the vanishing of the off-diagonal elements of this reduced density matrix at a 
succession of times is not identical with the decoherence of the corresponding histories, as 
will be discussed in Section IV. Yet the two notions are not unconnected. In the accessible, 
although unrealistic, model coarse grainings of the kind studied in this paper, where both 
notions are defined, typical mechanisms of decoherence ensure the validity of both. (See, 
e.g. [3], Section II. 6. 4) A general notion of mechanism of decoherence can be defined [19] 
that generalizes the reduced density matrix definition of decoherence in the context of the 
decoherence of histories and characterizes more precisely how they are connected. 
It would be clearer to use the terms "decoherence of histories" when referring to one 
notion and "decoherence of density matrices" when referring to the other. In this paper 
by "decoherence" we always mean the decoherence of sets of histories as defined precisely 
in this section. 



constitute the useful predictions in quantum mechanics. 

Among the coarse-grained decohering sets of alternative histories of this universe must 
be the sets that describe the quasiclassical domain of familiar experience. These consist 
of histories that, for the most part, are defined by ranges of values of "quasiclassical 
operators" correlated in time with high probability by classical phenomenological laws. 
We have discussed elsewhere the problem of distinguishing such quasiclassical domains 
from all other decohering sets of histories and in particular the problem of deriving the 
form of the "quasiclassical operators" that characterize them [2,33]. Here, we want to 
focus on a more specific and less general question. That is the question of the derivation 
and form of the classical phenomenological equations of motion assuming coarse graining 
is given that leads to decoherence. 

A simple form of coarse graining consists of averaging over some variables and following 
the remaining ones. Let us call the two classes of variables "ignored" and "distinguished" 
respectively. Widely occurring mechanisms of decoherence involve the rapid dispersal of 
quantum-mechanical phase information among the ignored variables as they interact with 
the distinguished ones. Thus, for example, sets of histories that distinguish the center- 
of-mass positions of bodies as light as a grain of dust can be very efficiently "decohered" 
by the collisions of the bodies with the omnipresent 3°K cosmic background radiation [7]. 
Such interactions can be expected to produce deviations from the predictability that char- 
acterizes classical behavior. That is, they produce noise. For quasiclassical behavior, such 
that the distinguished variables mostly resist the noise and follow approximately classical 
predictable paths, a high inertia is required for the distinguished variables. In general, 
therefore, we expect that a much coarser graining is necessary to achieve long stretches 
of predictable behavior than is needed to achieve mere decoherence, and decoherence, as 
we remarked, requires much coarser graining than is needed for mere consistency with the 
uncertainty principle. Furthermore, mechanisms that produce decoherence naturally lead 
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to processes, such as dissipation, that are necessarily included in the equations of motion 
that describe predictable behavior. There are thus connections among decoherence, noise, 
dissipation, and the amount of coarse graining necessary to achieve classical predictability. 
This paper explores those connections. 

The habitually decohering quasiclassical operators that characterize our everyday qua- 
siclassical domain include such "hydro dynamic" variables as the averages, over suitable 
volumes, of densities of energy, momentum, and other conserved or approximately con- 
served quantities. Such coarse grainings are not of the simple type described above in 
which the coordinates of a configuration space are separated once and for all into a set 
"distinguished" by the coarse graining and a set that is "ignored". First, the coarse grain- 
ings corresponding to these averages are not defined by ranges of coordinates. In addition, 
realistic coarse grainings are in general branch-dependent, meaning in this case, that the 
volumes over which the averages can be usefully taken are contingent on prior events in 
specific histories [2]. However, since the simple types of coarse graining are more easily 
analysed than the realistic ones, we shall begin our discussion with a class of model prob- 
lems that are based on distinguished and ignored coordinates and later return to how to 
generalize our results to more realistic situations. 

Our central result is a derivation of the classical equation of motion, including effective 
forces and noise, for a specific type of coarse graining of the histories of a class of quantum 
systems. Each system is assumed to have a Lagrangian that is the difference between a 
kinetic energy quadratic in the velocities and a potential energy independent of velocities 
but allowed to be fully non-linear. Coarse grainings are considered that distinguish a 
fixed subset of the coordinates of the system's configuration space while ignoring the rest. 
The initial density matrix of the closed system is assumed to factor into a product of a 
density matrix in the distinguished variables and another density matrix in the ignored 
variables. We show that when such sets of coarse-grained histories decohere, the quantum- 
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mechanical probabilities of the individual coarse-grained histories can be represented as the 
probabilities of the histories of a classical system evolving from distributed initial conditions 
under the action of a stochastic force. The initial conditions follow a Wigner distribution 
derivable from the initial density matrix. The distribution of the "total force" (including 
the inertial term "-ma") is a kind of generalization to histories of the Wigner distribution; 
it is derivable from the decoherence functional that measures quantum coherence. Neither 
of these distributions is in general positive. The effective classical equation of motion is 
the condition that the expected value of the total force vanish. When the noise arising 
from the stochastic force is almost negligible, so that the deviations from the effective 
classical equation of motion are small, then we achieve classical predictability. The noise, 
whether small or not, can be treated by incorporating a generalized Langevin force into 
the effective equation of motion. The distribution of the "total force" can also be regarded 
as a distribution of this Langevin force; it is, in general, non-Gaussian and, as mentioned 
above, not necessarily positive, as is to be expected from the quantum-mechanical nature 
of the problem. 

Our characterization of the effective classical equation of motion and the Langevin 
noise distribution leads to a method of identifying the equation of motion and a system- 
atic expansion procedure for calculating the noise distribution. These techniques are not 
restricted to linear problems or Gaussian noise. The key to the method is that the deco- 
herence functional, which depends on a pair of coarse-grained histories, can be expressed 
as a path integral over a quantity exp(L4), where we can expand A in a power series in 
the difference between the distinguished coordinates in one history and those in the other. 
For the familiar mechanisms of decoherence, to which we alluded above, the higher order 
terms in the expansion are expected to be small, and we can start by retaining only the 
linear and quadratic terms, which permit us to treat nonlinear equations of motion with 
Gaussian noise. The higher order terms then give the non-Gaussian corrections to the 



distribution of the Langevin force. 

As mentioned above, we would like to treat still more general and more realistic prob- 
lems, in which we escape the limitation to coarse graining that begins with distinguished 
and ignored coordinate variables, as well as the restriction to factored initial density matri- 
ces and the prohibition of velocity-dependent potentials. We discuss at length some ideas 
of how to free ourselves from these limitations. 

The details of the models are described in Section III. The various types of decoherence 
are discussed in Section IV. Section V introduces distribution functionals for the "total 
force" described above and the representation of quantum- mechanical probabilities in terms 
of them. In Section VI the equation of motion and distribution of noise are calculated for 
the well-known case of linear systems. The explicit generalization to non-linear cases is 
given in Section VII. The comparison of these results with the corresponding classical 
analyses is discussed in Section VIII. Section IX describes routes to more general coarse 
grainings and Section X contains some brief conclusions. 

There have been, of course, a great many discussions of the derivation of classical be- 
havior from quantum mechanics, and it is perhaps appropriate to offer a few comments on 
the similarities and differences between the present discussion and that great body of lit- 
erature. As already mentioned, we aim beyond the elementary discussions of the classical 
behavior of measured systems based on Ehrenfest's theorem, the WKB approximation, or 
the Wigner distribution. Such analyses do not usually treat noise, cover the effective classi- 
cal equations of motion including such phenomena as dissipation, consider coarse grainings 
(other than very obvious ones) or deal seriously with the probabilities of time histories by 
which roughly predictable quasiclassical behavior is inevitably defined. (Analyses based 
on the steepest-descent approximation to Feynman's path integral do consider histories 
but do not typically address the other issues.) 

Derivations of the equations of hydrodynamics from statistical physics, as, for example, 
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in [8,9,10, 41], necessarily include phenomena like dissipation. However, those accounts 
known to us derive equations of motion for the expectation values of hydrodynamic vari- 
ables. In quantum mechanics a system may be said to obey a classical equation of motion 
when the probability is high for the correlations in time that the equation of motion re- 
quires. For example, the center of mass of the earth can be said to obey Newton's law of 
motion when the probability is high that successive determinations of the position of the 
center of mass of the earth will be correlated according to that law. A complete derivation 
of classical equations of motion for quantum systems therefore requires the consideration 
of the probabilities for time histories, not just the study of the evolution of expected 
values. That is important, because the requirements of decoherence include the restric- 
tion on which sets of alternative histories may be assigned probabilities, while there is no 
such restriction on which expected values may be studied. Furthermore, it is through the 
study of the probabilities for histories that the probabilities for the inevitable deviations 
from classical predictability are most directly assessed. That is, only through a study of 
the probabilities of time histories can we accurately characterize the mixed quantum and 
classical-statistical noise that characterizes those deviations. 

The derivation of classical-statistical equations incorporating both classical determin- 
ism and stochastic noise has been extensively discussed for linear systems. Evolution 
equations for probability distributions on phase space were derived (in certain approxima- 
tions) from the evolution equation for the Wigner distribution by Caldeira and Leggett [11] 
and more recently by Hu, Paz, and Zhang [34] and by Zurek [35]. Langevin equations have 
also been extensively discussed for linear systems. Caldeira and Leggett [11], for example, 
discuss such equations, using techniques developed by Feynman and Vernon [12], and re- 
view earlier efforts. The extensive investigation of linear systems in the quantum optics 
literature is reviewed in [13] and treatments from the point of view of statistical mechanics 
can be found in [14, 41]. None of that work, however, explicitly considers the probabilities 
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of histories, by which classical behavior is necessarily defined. Neither is there explicit 
consideration of the decoherence of histories, which is a prerequisite for the calculation of 
probabilities in the quantum mechanics of a closed system. By using decoherence and the 
probabilities of histories we not only recover the standard Langevin equations for linear 
systems but also we can generalize these results to non-linear cases. 

III. Model Systems and Model Coarse Graining 

This paper is concerned with classical predictibility in a certain class of model quantum- 
mechanical systems described by particular classes of model coarse grainings. In this 
Section we shall specify these models and coarse grainings and review the theoretical 
framework in which the decoherence and the probabilities of histories are defined. We 
shall give a simplified version of this theoretical framework and later discuss how it is 
related to more general ideas in Section IV. That way the reader only interested in the 
models and not in these general connections can proceed immediately from Section III to 
Section V. We shall be brief. For greater detail the reader may consult [2, 33], and the 
references to the earlier literature found therein. 

We consider the quantum mechanics of a closed system in the approximation that gross 
fluctuations in the geometry of spacetime may be neglected. A background spacetime is 
thus fixed and, in particular, there is a well defined time, t. The usual apparatus of Hilbert 
space, states, operators, etc. then may be applied in the quantum mechanical process of 
prediction. The fundamental dynamics of the system is governed by an action S or its 
equivalent Hamiltonian H. The initial condition is specified by a density matrix p. 

The most refined possible description of a closed system makes use of a set of fine- 
grained histories. The most familiar set of fine-grained histories are the possible paths, 



Recently, Dowker and Halliwell [40] have studied the decoherence of histories in explicit 
linear models and also, in effect, derived their classical equations of motion, although not 
including a description of dissipation and noise. 
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q ' {t), in a configuration space of generalized coordinates that completely describe the 
system. These paths are the single-vaued functions q^(t) on a fixed time interval, say 
[0,T]. For example, in a system of scalar fields each q^(t) might be the value of a field 
at a different spatial point considered as a function of time. (Our description is thus 
not restricted to non-relativistic physics.) Configuration space fine-grained histories are 
the starting point for sum-over-histories formulations of quantum mechanics and for the 
model coarse grainings we shall mostly consider in this paper. More general possibilities 
are discussed in the next Section. 

A partition of an entire set of fine-grained histories into exhaustive and exclusive classes 
defines a set of coarse-grained histories; each class is an individual coarse-grained history. 
The individual coarse-grained histories in an exhaustive set may be grouped into new 
exclusive sets. That is an operation of further coarse graining, yielding a coarser-grained set 
of alternative histories. The inverse operation of dividing a set of coarse-grained histories 
into smaller classes of the fine-grained histories is an operation called fine graining. Sets 
of coarse-grained histories are partially ordered with respect to the operations of fine and 
coarse graining. A rich variety of coarse grainings is possible. As we mentioned before, the 
histories of the quasiclassical domain of everyday experience, for example, are defined by 
coarse grainings utilizing ranges of values of averages over suitable spatial regions of such 
"hydro dynamic" variables as the densities of energy, momentum, charges, and currents. 
Such realistic coarse grainings are in general branch-dependent, that is, contingent on 
which of many possible events have happened. 

Because of their branch dependence, and because of their indirect relation to the fun- 
damental fields, realistic coarse grainings are not as theoretically tractable as some model 
coarse grainings that can be studied. In this paper we shall study a familiar and instruc- 
tive class of model coarse grainings in which the coordinates, q ' : of configuration space are 
divided into ones x a that are distinguished by the coarse graining and the remaining ones 
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Q , which are ignored. (For example, in a simplified model of a universe of particles, the 
x a might label the center-of-mass positions and orientations of a group of massive bodies 
such as the planets and the Q would then be all the rest of the coordinates, including 
the internal coordinates of the bodies' constituents and the coordinates of gas molecules, 
etc. that interact with the planets.) Coarse-grained histories of this type are labeled by 
partitions of the paths x a (t). An individual coarse-grained history consists of a path x a (t) 
along with all possible paths Q (t). 

Further coarse graining of the classes of histories labeled by the paths x a (t) can be 
defined by sets of intervals exhausting the whole range of the x a at a discrete sequence 
of times t\ < ti < ■ ■ ■ < t n . These correspond to intervals on the entire configuration 
space of q" that are unrestricted in the Q but are subject to the stated restrictions on 
the x a . We denote the exhaustive sets of such intervals at the successive times by {A3, }, 
{A^, 2 }, • • •, {A" n }. The index k on {A*} labels the particular set, a^ labels the particular 
interval in the set, and t^ is the time. An individual coarse-grained history in such a set 
consists of the paths, q (£), that thread a particular sequence of regions, e.g. A 5, at t\, 
A^, 2 at £2, etc, and the whole set of coarse-grained histories is exhausted as the different 
possible ways which paths may pass through the regions are enumerated. A particular 
coarse-grained history thus corresponds to a sequence (cq, • • • , a n ), which we shall often 
abbreviate as just a. 

The coherence between individual histories in a coarse-grained set is measured by the 
decoherence functional. This is a complex functional defined on all pairs of coarse-grained 
histories in an exhaustive set of such histories. The sets of histories under discussion are 
coarse grainings of fine-grained histories that are paths q' (t) on the time interval [0, T]. 
The sum-over-histories formulation of quantum mechanics is, therefore, convenient for 
introducing the decoherence functional; more general formulations are discussed in the 
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next Section. For the coarse grainings under discussion, the decoherence functional is 

D(a', a) = ISq' Sq5(q' f - q f ) expji (%'(r)] ~ %( r )]J / n }p{%^ Qo) ■ ( one ) 

The path integral over q (t) is over all paths that start at q^ at t = 0, pass through the 
intervals A^, A^ 2 , • • • , A" n of the x a at times t\ < • • • < t n , and wind up at ql at time 
T. The path integral over q'^(t) is similarly defined. The integrals include an integral 
over the initial q% and final ql. Here, piq^qo) is the initial density matrix in the q 
representation and ^[(/(r)] is the fundamental action. The measure for the path integrals 
is the standard one induced by the Liouville measure in phase space. It is described 
explicitly in Appendix (a). Equation (one) has been compressed by omitting the indices 
on the q ' (t) and by denoting the entire sequence (cq, • • • , a n ) at (t\, ■ ■ ■ , t n ) by a single 
index a. We shall employ similar conventions in the rest of the paper. 

The coarse grainings under consideration distinguish only the x a in the division 
q = (x a , Q ). Following an analysis of Feynman and Vernon [12], the integrals over the 
Q may therefore be carried out over their whole ranges unrestricted by the particular 
coarse-grained histories considered. Suppose that the action may be decomposed as 

S[q(r)} = S ircc [x(r)} + S [Q(r)} + S int [x(r), Q{r)\ . (two) 

(The use of the subscript "free" does not mean that there is no potential energy term 
in <Sf ree . There is, in general. Rather it means that the action of the x's is free of any 
interaction with the Q's.) The integral over the Q's defines W - - a functional of the paths 
x'(t) and x(t) and a function of their initial endpoints Xq and xq -- as follows: 



exp (iW[x'(T),x(T);xQ,xo)/h) p(x'q,xo) = SQ' 6Q5(Q'f-Qf) 
x exp\i(So[Q'(r)] + S [nt [x'(r), Q'(t)} - S [Q(r)} - S- mt [x(r), Q(r)])/h 



xp(x'q,Qq;x ,Q ) . (three) 
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Here we have introduced the reduced density matrix p = Sp p associated with the coarse 
graining, 

p(x' Q , xq) = / dQo p(x' , Qo; ^0, Qo) • (four) 

The functional W is only a slight generalization of the Feynman- Vernon influence phase 
and we shall continue to call it that. It depends on the endpoints x'q and xq implicitly 
through the paths x (t) and x(t) because the actions on the right hand side of (three) 
are functionals of these paths. There is also an explicit dependence on Xq and xq arising 
from the dependence of p on these variables in (three) and (four). We use the notation 
W[x'(t),x(t);x'q,xq) to indicate this dependence, the square bracket to indicate depen- 
dence on functions and the round bracket to indicate dependence on variables, and we 
maintain this notation for other cases. The quantity W depends on the time interval T as 
well, but we have not indicated this explicitly. 

The decoherence functional may then be expressed in terms of W and p as follows: 

D(a',a) = I Sx' / 8x 8(x'f — xA 

Ja' Ja 
xexpji^free^^r)] - 5free[x(r)] + W[x'(r), x(r); x' , x )J /hjp(x , x ) . (five) 

Thus all the contribution from the ignored variables is summarized by the functional 

W[x'(t), x(t); x'q, xq). 

Restrictions on the form of the actions S^ee t^] ; Sb[Q]> an< ^ ^intt^? Q]i as weu as on the 
form of the density matrix p, will be needed for the explicit derivation of the equation of 
motion, as described in the subsequent sections. We shall, for example, generally make 
the usual assumption that Sf Tee [x(r)] has a simple "kinetic minus potential" form 

1 



T 

(six) 



-x\t)Mx(t)-V(x(t)) 



5free[a;(r)] = dt 
JO 

where we have used an obvious matrix notation x'Mx = T, a i ) x a M a i ) x . Similar assumptions 



will be made for Sq[Q(t)]. We shall assume that Si n i[x(r), Q(t)] is local in time, that is, 
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of the form 

Swt[x(r),Q(r)]= dt L int [x(t),Q(t)} , (seven) 

JO 

with Li n t independent of velocities. These assumptions restrict us only to a widely ap- 
plicable class of models and it is likely that similar results can be obtained from weaker 
assumptions. 

A more restrictive assumption concerns the form of the initial p. We shall make the 
conventional assumption [12, 11] that p factors into a density matrix in x and another in 

Q 

p(x' Q , Q' ; xq, Qq) = p(x' Q , x Q )p B (Q' Q , Qq) , (eight) 

so that the variables distinguished by the coarse graining are initially uncorrelated with 
those it ignores. Of course, we do not necessarily expect initial density matrix of the 
whole universe to factor as in (eight), but, for widespread mechanisms of decoherence 
that operate essentially locally in space and time when compared with cosmological scales, 
(eight) is an excellent approximation. For example, scattering by the cosmic background 
radiation can efficiently decohere alternative positions of the center of mass of a massive 
body coarse grained on centimeter scales [7] . The coordinates of the body and radiation 
may be correlated in the wave function of the universe, but on the local scales where the 
mechanism operates they are effectively uncorrelated, as described by (eight). 

Factorization has a number of helpful consequences for the form of W and p defined 
in (three) and (four). First, p is given by 

p(x' 0: x ) = p(x' ,x ) . (nine) 

Second, and most importantly, the influence phase W contains no explicit dependence on 
x'q and xq and we may write W[x (r), x(t)]. This will simplify the form of the equation of 
motion, which would otherwise contain terms arising from the explicit dependence on x'q 
and xq. 
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For most of this paper, therefore, we are considering a class of models defined by 
coarse grainings that distinguish a fixed subset of the coordinates of configuration space, 
by actions that have the simple forms (two), (six), (seven), and by an initial density 
matrix that factors as in (eight) . Especially simple examples of such models are the linear 
oscillator models studied by Feynman and Vernon [12] and by Caldeira and Leggett [11]. 
In these models a distinguished oscillator is coupled linearly to a large number of other 
oscillators constituting a thermal bath characterized by a temperature Tg. The density 
matrix is assumed to factor as in (eight). Let x(t) be the coordinate of the distinguished 
oscillator and u)r its frequency renormalized by its interactions with the others. Then, in 
the case of a continuum of oscillators, cut off at frequency O, and in the Fokker-Planck 



limit of /cTg >> hQ » hujR, Caldeira and Leggett find for the influence phase 

7/ 

'0 h jo 
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W[x (r), x(t)] = —M'j / dt [x x — xx + x'x — xx~\ +i / dt [x (t) — x(t)] 



(ten) 
were 7 is a coupling constant summarizing the interaction of the distinguished oscillator 
with the rest. 

IV. Decoherence 

(a) Decoherence in General 

Quantum mechanics predicts the probabilities for the individual members of a set 
of alternative coarse-grained histories only when there is negligible quantum-mechanical 
interference between the individual members of the set [5, 6, 2]. Only then do the squares 
of amplitudes define probabilities that are consistent with the sum rules of probability 
theory. Sets of histories that exhibit negligible interference are said to decohere weakly 
[33]. 

However, in quantum mechanics we are not interested just in sets of histories that 
are consistent in the sense that they can be assigned probabilities satisfying probability 
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sum rules. We are interested also in sets of histories that constitute the quasiclassical 
domain of everyday experience. It is this quasiclassical domain that lies at the root of 
the interpretation of quantum mechanics. It is through an understanding of this domain 
that quantum mechanics acquires utility for our experience; "measurement" situations arise 
precisely when variables become highly correlated with the quasiclassical domain. Stronger 
notions of decoherence are therefore useful to characterize the realistic mechanisms of 
decoherence that lead to a quasiclassical domain [33]. In this section we shall review 
several notions of decoherence that we have described in previous work, and we shall 
discuss their connections with each other and with some other notions of decoherence that 
have been introduced in the literature. We shall be brief. For greater detail the reader 
may consult references [2] and [33]. 

We begin by recalling how a set of alternative, coarse-grained histories of a closed 
system is described in quantum mechanics. The simplest kinds of histories are specified 
by giving independent sets of alternatives at a sequence of times t\ < t<i < ■ • • < t n . 
In the Heisenberg picture, alternatives at one moment of time t^ correspond to a set of 
projection operators {Pa k (tk)}. The index k denotes the set of alternatives at time t^, 
while the index a^ denotes the particular alternative within that set. These projections 
represent exclusive alternatives, so they are orthogonal for different alternatives, and they 
represent an exhaustive set, so they sum to unity over all alternatives. For example, for the 
coarse graining by ranges of a distinguished set of variables x a described in the preceding 
Section, the projection Pj% (t^) would just be the projection onto the range A^ fe at time t^. 
An individual history in a set defined by a sequence of such sets of alternatives corresponds 
to a sequence of particular alternatives a = (cq, • • • , a n ). Each history is represented by 
the corresponding chain of projection operators 

C a = PS n (t n ) ■ ■ ■ P^) . (fourone) 

A completely fine-grained set of histories would consist of one-dimensional projections 
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onto complete sets of states at each and every time. Sets of histories defined by sets of 
projections that are not all one-dimensional or not at every moment of time are said to be 
coarse-grained. 

In the Heisenberg picture, every exhaustive set of orthogonal projection operators in 
Hilbert space {Pa k } represents, at any time, some set of alternatives for the system. Of 
course, the alternatives corresponding to a given set of projection operators will have differ- 
ent descriptions in terms of fundamental fields when different times are assigned to them. 
Similarly a sequence of sets of projection operators define different alternative histories 
when different times are assigned to the sets of projections. However, an assignment of 
times leads to meaningful alternative histories only if the ordering of the times corresponds 
to ordering of the projections as in (fourone). If this time ordering is not respected, two 
inconsistent sets of alternatives could be assigned the same time and the resulting alter- 
natives would not be meaningful. 

While in quantum mechanics we usually consider sets of histories consisting of inde- 
pendent sets of alternatives at sequences of times, a more realistic description is achieved 
by generalizing this notion in two related ways [2, 3]. First, if a = (a\, ■ ■ ■ , a n ) is a his- 
tory and j3 = (/3i, • • • , (3 n ) is a distinct history then we may consider the coarser-grained 
alternative that the system followed either history a or history (3. The alternative a or (3 
is represented by the sum of the chains for a and (3 

C a or (3 = C a + Cp (fourtwo) 

but is not itself necessarily a chain of projections of the form (fourone). More precisely, 
if {a} is a set of alternative histories for the closed system defined by sets of alternative 
projections at a sequence of times, then coarser-grained sets are defined by partitions of 
the {a} into exclusive classes {a}. The classes correspond to the individual histories in 
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these coarser-grained sets and are represented by operators C a that are sums of the C a 

Ca = Y._C<*= E P S n (t n ) ■ ■ ■ P^(tl) (fourthree) 

aea (ai,-,a„)ea 

These generalized C a need not themselves be chains of projections and thus we sometimes 
extend the use of C a to denote the operators representing individual histories in these 
more general coarse-grained sets. 

The second important generalization is to allow the histories to be branch dependent, 
that is, for the set of alternatives at time t^ to depend on the values of earlier labels 
cq, • • • , CKfc_i [2, 6]. Branch dependence is important, because in a quasiclassical domain 
past events may determine what is a suitable quasiclassical variable. For example, if a 
quantum fluctuation gets amplified so that it leads to condensation of a galaxy in one 
branch and no such condensation in other branches, then the outcome clearly influences 
what are suitable quasiclassical variables in the region where the galaxy would form. 

For the general case of branch dependence, a better notation than (fourone) for chains 
of projections would be the following: 

C a = Pa n (t n ; a n -l, •••,«!) P u~\ (*n-i; a n-2, • • • , «l) • • ■ P^ih) , (four four) 

where the Pg, (t^; a^_i, • • • , cq) define an exhaustive set of orthogonal projection operators 
as «£, runs over all values for fixed cq, • • • , a^-i, corresponding to an exhaustive set of 
mutually exclusive alternatives for the closed system. In limiting ourselves to projections 
that depend only on previous alternatives rather than future ones we have incorporated a 



In Ref. [2], Section X, we unnecessarily eliminated the possibility of branch-dependent 
chains of the form (fourfour) and restricted attention to sets of ct^'s that were indepen- 
dent of one another. We did that in order to safeguard a special derivation of the weak 
decoherence condition for assigning probabilities to alternative coarse-grained histories. In 
fact, that special derivation is unnecessary, and as we shall see below weak decoherence 
can easily be seen to be the necessary and sufficient condition for the probability calculus 
to apply to histories that are chains or sums of chains, whether or not the choice of the 



set P^ k is branch-dependent. 
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notion of causality consistent with the usual arrow of time in quantum mechanics. Further 
generalizations to formulations without an arrow of time are possible [24] . 

In the models that we treat in detail in this article, we do not make use of branch de- 
pendence, since we assume a fixed division of coordinates into those distinguished, x a , and 
those ignored, Q , and the A's (intervals of x values) are taken to be branch-independent. 
We can thus employ the simplified notation (fourone). It should be borne in mind, though, 
that realistic coarse-grained histories do involve branch dependence. 

Having in hand this discussion of the possible sets of alternative coarse-grained histories 
of a closed system, we can now turn to the various notions of their decoherence. The central 
quantity, the decoherence functional, is defined generally for pairs of histories in a coarse- 
grained set by 

D(a, a) = Tr[C a ipCl] (four five) 

for a density matrix p representing the initial condition and operators C a representing the 
individual histories. 

The necessary and sufficient condition for probability sum rules to be satisfied is 

ReD(a f , a) = , a ^ a . (foursix) 

In previous work [33] we have called this the weak decoherence condition to distinguish 
it from stronger decoherence conditions we shall discuss below. When a set of histories 



As pointed out to us by Bob Griffiths, we have incorrectly attributed the weak decoherence 
condition (4.6) to him (in reference [2] and elsewhere). In fact, Griffiths [15] and Omnes 
[6] employ a weaker condition than (4.6) as the necessary and sufficient condition for the 
"consistency of histories". That is because they require, in the notation of (fourone), that 
two chains of projections C a and C a > must have the real part of their interference term 
vanish only if C a + C a i is another chain of projections, whereas we require it in the case of 
all the chains C a and C a i (see below). Our weak decoherence applies to coarse grainings 
that are allowed to unite any two histories in the set being studied, while the "consistent 
histories" condition of Griffiths and Omnes applies only to some of those coarse grainings. 
For example, if a and a differ in only one index, our approach and theirs give the same 
conditions, but if a and a' differ in more than one index, then Griffiths and Omnes do 

not always require that Re^r^^pC^)} = but rather a weaker condition. In essence 
Griffiths and Omnes restrict themselves to histories defined by independent alternatives 
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weakly decoheres, the probability of a history a is the corresponding "diagonal" element 
of the decoherence functional 

p(a) = D(a, a) . {four seven) 

Eq. (fourone) is the necessary and sufficient condition that these numbers obey the sum 
rules of probability theory. All that is needed to show this is to notice that the probability 
that either of two histories will happen is the sum of the probabilities of the two individual 
histories if and only if the sum of the interference terms represented by (foursix) vanishes. 
Weak decoherence is the criterion by which quantum mechanics discriminates between 
those sets of histories that can be assigned probabilities and those that cannot. 

A stronger notion of decoherence is provided by the medium decoherence condition [33] 

D(a , a) = , a ^ a . (foureight) 

Clearly medium decoherence implies weak decoherence, but not the other way around. 
Medium decoherence is a consequence of realistic mechanisms that are widespread in the 
universe. It is, therefore, a natural condition to impose in characterizing a quasiclassical 
domain. It is not, by itself, sufficient to single out a quasiclassical domain. The necessary 
further criteria are a large part of the subject of this paper. 

When the initial state is pure, exact medium decoherence is equivalent to the existence 
of generalized records for each history in the decohering set [33] . To see this, notice that 
for a pure initial condition, p = |\1/)(\1/|, 

D(a',a) = ((#|CJ) (C a /|*)) . (fournine) 



at a sequence of times. Each history corresponds to a sequence of such alternatives and 
therefore to a chain of projections. They do not therefore incorporate branch dependence, 
at least in the sense of (fourfour). In the models studied in this paper, the stronger 
conditions are, in fact, satisfied. The whole of any off-diagonal element of the decoherence 
functional approximately vanishes - - not just the real part. Also, when more of the 
indices in a and a! differ, the decoherence condition (4.8) is satisfied more strongly, not 
more weakly. This gives us some confidence that our stronger conditions are physically 
realistic. 
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To every set of alternative histories there corresponds a resolution of the pure initial 
state into branches 

|*> = £( C al*» = £„ n P% n (t n ) ■ ■ ■ Pi(ti)|*> • (fourten) 

a 

For a set of histories obeying exact medium decoherence, the branches are orthogonal, 
as (foureight) and (fournine) show. Therefore, there exists at least one set of orthogonal 
projection operators {R a } that project onto these branches 

R a \^) =C a \®) , (four eleven) 

where 

R a R a ' = fiaa'Ra and y^R a = 1 , (fourtwelve) 

a 

so that the projections {R a } are exclusive and exhaustive. 

Now we can also assign projection operators R a to coarse grainings of the histories 
{a}, that is, sums of chains of projections {C a }. To every such coarser-grained history, 
representing the union of a subset of the histories {a} or the sum of the corresponding 
{Ca\, we assign the projection operator that is the sum of the corresponding {R a }. That 
is perfectly consistent with the extension of (foureight) to the coarser-grained histories. 
The resulting K's have the property that progressive fine graining of the coarser-grained 
histories results in a sequence of nested i?'s, projecting onto smaller and smaller subspaces 
of Hilbert space, where each such space is a subspace of the preceding one. 

When the branches C a \^>) 7^ do not form a complete set of orthogonal states for 
the Hilbert space, there can be many sets of projections {R a } that obey the conditions 
(foureleven) and (fourtwelve). When the branches do form a complete set, then the {R a } 
are unique; they are just the projections onto the single states C a \^). The set of histories 
{a} is then said to be full [33]. The i?'s for coarse grainings of the histories {a} are then 
also unique: for the further coarse-grained history that corresponds to the union of a given 
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subset of the {a}, the corresponding R is just the sum of the relevant R a and projects 
onto the space spanned by the corresponding vectors C a \^>) = R a \^). 

When assigned a time after t n in the sequence t\ < t < • • ■ < t n , the -R a 's may be 
thought of as representing generalized records of the histories. They may not represent 
records in the usual sense of being constructed from quasiclassical variables accessible to 
us, but the condition (foureleven) means that at any time there is complete information 
somewhere in the universe about the histories {a}. 

The above discussion shows that medium decoherence in a pure initial state implies 
the existence of generalized records. The converse is also true. The existence of orthogo- 
nal generalized record projections satisfying (foureleven) ensures the medium decoherence 
of the corresponding set of histories through (foureight) and (fournine). Exact medium 
decoherence can thus be characterized by records, and the physical formation of records is 
a way to understand mechanisms by which medium decoherence occurs. In the example 
implicit in the work of Joos and Zeh [7] , histories describing successive alternative posi- 
tions of a dust grain, initially in a superposition of positions about a millimeter apart, very 
accurately satisfy the condition of medium decoherence simply by virtue of the scattering 
by cosmic background radiation photons. The successive scatterings of these photons ef- 
fectively create records of the histories of positions of the dust grain in the electromagnetic 
degrees of freedom. The commuting records of successive positions are stored indepen- 
dently in the vastness of cosmological space as the photons move off at the speed of light. 
They may not be accessible to us, but their existence is a way of understanding how this 
mechanism of medium decoherence works. 

The permanence of the past is a feature of the quasiclassical domain that is naturally 
explained by medium decoherence when there is a pure initial condition. By permanence 
of the past we mean the feature of a quasiclassical domain that what has happened in 



For a more detailed analysis in terms of histories, see [3]. 
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the past is independent of any information expressed by a future projection. Neither 
the decoherence of past alternatives nor the selection of a particular past alternative is 
threatened by new information. "The Moving Finger writes; and having writ, Moves on: 
nor all thy Piety nor Wit Shall lure it back to cancel half a Line, Nor all thy Tears wash 
out a Word of it" [38]. 

In other words, we are discussing the property of the decohering coarse-grained histories 
{a} that, at any of the times t k , there is, for each history up to that time, an effective 
density matrix [3] 

Pg k (tk)---Pi(t 1 )pP^ 1 (t l )---Pl(t k ) fourthirteen) 

Tr[pk k (t k )---Pl 1 (t 1 )pPl 1 (t 1 )---Pk k (t k )} ' U ° Ur %U<< " ] 

which can be utilized for all future predictions without concern about the other values of 
«1, • • • , otfc. For a pure state, this corresponds to the "reduction of the state vector" 

|*) ->W-2P* fc (**)•• -Pi(*l)|*) , (four fourteen) 

where N is the trace in (fourthirteen). 

When the past is permanent, we may still lose the ability to retrodict the probabilities 
of alternatives in the past through the impermanence or inaccuracy of present records but 
not from the failure of those past alternatives to decohere in the face of the projections that 
describe information we acquire as we advance into the future. Yet we know that such 
continued decoherence of the past is not guaranteed in general by quantum mechanics. 
Adjoining future alternatives to a set of histories is a fine graining of that set and in 
general a fine graining of a decoherent set of histories may no longer decohere. Verifying 
the continued decoherence of all the past alternatives as we fine-grain our set of histories 
to deal with the future would in general require significant computation. We would have to 
check that the branches corresponding to every alternative past that might have happened 



For further and less informal discussion see [3], Section II. 3. 2. 
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continue to be orthogonal in the presence of their newly adjoined sets of projections. Yet 
we adjoin sets of projections onto ranges of quasiclassical operators without making this 
calculation, secure in the faith that previous alternatives will continue to decohere despite 
this fine graining. It is this assumption of continued decoherence of the past that permits 
the focus for future predictions on the one branch corresponding to our particular history 
and the discarding of all others. In other words, we pointed out above, it is the permanence 
of the past that permits the "reduction of the state vector" . 

If we consider, instead of the set of histories {a} = {(cti, • • • , a n )}, the set of abbrevi- 
ated histories {(a\, • • • , a^)} with k < n, running up to time t^ < t n , that is an example of 
a coarser graining of the set of coarse-grained histories {a}. As these abbreviated histories 
get further a's (a^+i then a^+2) etc.) adjoined to them, that represents a sequence of fine 
grainings of the coarser- grained histories {(cq, • • • , a^)}. 

If p is pure, there are nested records R corresponding to these abbreviated histories. 
When further cc's are adjoined, as the histories unfold from (cq, • • • , a^) to (a\, • • • , (*k+l)-> 
etc., up to (a\, • • ■ ,a n ), the nested record operators are projections onto subspaces of 
Hilbert space that progresively narrow. In such a situation, the past always continues 
to decohere as the histories advance into the future. In fact, the physical formation of 
nested, generalized record operators R guarantees, in general, the permanence of the past, 
including not only the permanence of its decoherence, but also the permanence of the 
selection of particular past alternatives as well. 

If we relax the condition of exact decoherence and consider approximate medium de- 
coherence, defined by the approximate validity of eq. (foureight), then it is possible to 
understand more about the formation of records and the origins of medium decoherence 
for the kind of model coarse grainings studied in this paper, which distinguish particular 
coordinates and ignore a large number of others. The Hilbert space is a tensor product of 
a Hilbert space of functions of the distinguished coordinates, 7i x , and one for the ignored 
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coordinates, 7v* . The coarse-grained histories consist of chains of projections at times 
t\, t2, • • • that, in the Schrodinger picture, act only on H x and represent partitions of some 
complete set of states in 7i x at each time. In the models, these states are, in fact, just the 
localized states in x a at each time, but for greater generality we shall consider different 
complete orthogonal sets at each time represented by wave functions {4>ri( x )} at time t\, 
{(pr 2 ( x )} & t time t^ etc. In the subsequent discussion we shall augment this notation to 
indicate the branch dependence of the possible orthogonal sets. 

In order to make a connection with the thinking of some authors, we shall take a brief 
excursion into the Schrodinger picture, while assuming a pure initial state represented by 
a wave function ty(x, Q, to). Let us follow for three steps the Schrodinger evolution of this 
initial \1/. It may be evolved to the time t\ of the first set of alternatives and expanded in 
the first complete set {(p}.(x)} as 

*(z, <2, to) -»• Y ^i ( x )* r i (Q) • (four fifteen) 

n 

The right-hand side is *&(x,Q,ti) in the Schrodinger picture and {x(t\) = x,Q(t{) = 

Q\ty) in the Heisenberg picture. The coefficients Xn are in general neither normalized 

nor orthogonal for different r\. A coarse graining divides the r\ up into exclusive and 

exhaustive sets {cci}. The sum on the right of (fourfifteen) may be similarly decomposed: 

V(x,Q,ti) = ^2^2 ^ri (^)Xri (<3) • (four sixteen) 

u\ r\€(X\ 

The result of the evolution of each branch J^ ^ ) r 1 ( ;r )Xri(Q) m (foursixteen) to the 

rieai 

time £2 °f the next set of alternatives may again be expanded in a new complete set of 
functions in the rr's. In general this set will be branch dependent, that is dependent on 
a\. We use the notation {4>r 2 1 ( x )} to indicate this, and the expansion then has the form 

Yl ^i( x )Xr 1 (<5) -»■ Yl Yl ^(^Xr^iQ) ■ (four seventeen) 

rieai a 2 r 2 ea 2 
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Similarly for the evolution from t<i to £3: 

E cj> 2 r ^(x) X ^(Q) - E E < a2ai (*)Xr 3 2ai (Q) > {four eighteen) 

T2ta2 ct3 r 3 ea 3 

and so forth. Thus, 

*(x,g,t 3 )= E ( E $T ai ( x )Xrt ai (Q)\ ■ {fournineteen) 

a 3 a 2 ai \r3ea3 / 

The term in the brackets is the branch wave function corresponding to the history 

a = (a\, 0.2, 0:3). The overlap of branches gives the decoherence functional according to 

(fournine) and (fourten). Thus 

^(«» = E E (^ 2ai aI'' 2 "' 1 ) {xT 1 > X^' 1 ) ■ (fourtwenty) 

r' 3 ea' 3 r 3£0- 3 

If the complete set at time £3 is independent of previous alternatives, so that it is not branch 
dependent, then the decoherence of the last alternative is automatic. The scalar product 
(0p,, <j> , ) = S r r i and so the right-hand side of (fourtwenty) is diagonal in r 3 . Otherwise it 
is a non-trivial condition because {^v" 2 " 1 } an d {4> r , 3 2 1 } are different orthogonal sets when 
a's do not coincide. Suppose the Hilbert space 7v* is effectively very large compared to 7i x , 
in the sense that at each time t^ a one-to-one correspondence is established between 7i x 
and a tiny portion of liH . Then we may expect the scalar products between x's differing 
by any index to be typically very small , leading to small values of D(a, a') for a' 7^ a. 
The summations over unorganized phases between different values of r 3 may make the 
off-diagonal elements of D(a , a) even smaller. That is approximate medium decoherence. 
In the approximation in which all the x's are really orthogonal to one another, the records 
R a that accomplish projections onto branch wave functions such as those in (fournineteen) 
may then be taken to be projections onto the corresponding sets of x's: {Xr 3 2ai | r 3 eQ; 3}- 



Of course, as more and more times are added to the histories, and the set of ct's grows 
larger accordingly, we expect that eventually even the space liP will be exhausted and x's 
no longer orthogonal. 



28 



The reduced density matrix in x may also be constructed. It is 

p&>*) = £ £ £ £ ^ a V)4 a2ai *(*) (x?r , 4 ar 

(four twenty one) 
involving a sum over all branches. If indeed all the x's are approximately orthogonal to one 
another because of the largeness of the Hilbert space Tir* , then p(x' ,x) is approximately 
diagonal in the histories {a}. 

Many authors have considered, not the full density matrix (fourtwentyone) , but the 
reduced effective density matrix for each branch [cf. (fourthirteen)], relevant when the 
decoherence of earlier alternatives may be assumed. The applicable portion of p(x , x) 
[(fourtwentyone)] is 

P a *°* y iX ) = ££ £ £ 4% 2ai M 4T ai *W {xT 1 > x% ai ) ■ 

"3 a' 3 r' 3 ea' 3 r 3 ea 3 

(fourtwentytwo) 
Approximate orthogonality of these x's for different 03 leads to approximate diagonality 
of p a2CXl in the ct3's. However, a much stronger condition on the x's is needed to ensure the 
decoherence of whole histories as in (fourtwenty) than is needed to ensure the diagonality 
of (fourtwentytwo). We need the approximate orthogonality of the {Xr^ 1 } n ot only when 
the r3 belong to different sets C13, but also between functions corresponding to different 
values of cq and «2- The strength of these conditions may be appreciated by noting two 
facts: First, since any coarse graining of a decoherent set is also decoherent, the x's must 
be approximately orthogonal in all previous steps. That is, the {x^ 1 } i n the example must 
be approximately orthogonal for different cq's and when the r<i lie in different a^'s- There 
is a similar condition for the {xri}- Since x and Q interact between the times t\, • • • , t n 
none of these conditions is a simple consequence of the others. Second, nothing in our 
discussion has fixed the choice of the times t\, • ■ • , t n . Of course, they must be sufficiently 
separated for the interactions to disperse the phases, but once that is satisfied, we expect 
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decoherence to hold for a range of times giving rise effectively to even more conditions. 
The decoherence of histories is a much stronger requirement than the diagonality of density 
matrices. 

Some authors [17, 42, 18] have discussed how, in the case of a pure initial density 
matrix ^>(x', Q')^*(x, Q), one can use the Schmidt decomposition to achieve some of the 
objectives of decoherence. In the notation we have been using for the Schrodinger pic- 
ture in eqs. (fourfifteen) - (four nineteen), we can choose, at each time t^, the functions 
4> rk k ~ 1 1 (x) to be Schmidt functions, which means that their coefficients Xr^ 1 1 (Q) 
are orthogonal for different values of r^ (although still not normalized). In that way the 
reduced density matrix for each branch at time t^ becomes diagonal in r^, as we can see 
for the case k = 3 in (fourtwentytwo). It is then, of course, also diagonal in a^ (0:3 for 
the case k — 3). However, this Schmidt construction does not guarantee the decoherence 
of even the final alternative as defined by the decoherence functional (fourtwenty) (unless 
we have a\ = a^ and «2 = a 2 ). That is because the Schmidt functions are necessarily 
branch- dependent and the orthogonality of the %'s for one branch does not guarantee the 
orthogonality between different branches that would lead to decoherence of the final alter- 
native. In any event, the main thing for the decoherence of histories is to have D(a , a) 
diagonal in a.<i an d ol\ as well. That is not guaranteed at all by the Schmidt procedure, 
which does not imply the orthogonality of \ functions for different values of the ct's. In 
particular, the Schmidt procedure does not guarantee the permanence of the past discussed 
above. The relation between diagonalization of the reduced density matrix at successive 
times and the decoherence of histories will be further discussed in [19]. 

Now let us consider the situation with an impure initial state represented by a density 
matrix p. An impure initial state could be fundamental, representing the initial condition 
of the universe. However, even if the cosmological initial condition is pure, a reduced 
density matrix, in which some coordinates already have been traced over, can be a use- 
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ful description of local physics. A familiar example is the cosmic background radiation. 
Imagine that the background photons have all been pair produced in a pure cosmological 
initial state. In that case, for every photon near us, there would be a correlated photon 
with equal but opposite momentum on the other side of the universe. The local physics, 
however, would be accurately described by a nearly thermal reduced density matrix in 
which the distant photons had been traced over. 

In [33] we discussed a "strong decoherence" condition that is the analog of (foureleven) 
for density matrices. A set of histories was said to decohere strongly when, for each history 
in the set, there exists record projections satisfying (fourtwelve) such that 

C a p = R a p . (four twenty two) 

However, we shall see that for highly impure states this is too strong a condition to char- 
acterize usefully a quasiclassical domain. 

We can express an impure density matrix p in terms of its eigenstates and eigenvalues: 
p = ^2 ^l*^) (*1 > (f ourtwentythreea) 

or, in the special case of x and Q variables, 

p (x', Q'; x, Q) = V -k^^(x' ^Q 1 )^*^ (x,Q) • (f our twenty threeb) 

A* 

Here the tt^ are the probabilities of the initial states |\I/^). 

If p has non-zero probabilities for many states, then (fourtwentytwo) is difficult to 
satisfy. If p has non-zero probabilities for an orthogonal set of states |\I/^), then (fourtwen- 
tytwo) would imply 

C a \^) = R a \^) for all \V») . (four twenty four) 

If the set |\1/^) is complete, then (fourtwentyfour) implies that C a = R a . That can be 
satisfied only in the trivial case in which all the P's for all times commute with one 
another. We shall henceforth ignore "strong decoherence" for impure density matrices. 
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We could define, when p is impure, a kind of "medium strong decoherence" , in which 
we would have, for each |\1/^) with non-zero probability 7r^, a generalized record projection 
operator R& such that 

C a \^) = i2g|tf") , (fourtwentyfive) 

where for each \i the R& are exclusive and exhaustive projections. This requirement would 
mean medium decoherence separately for each |\1/^) (with tt/j, ^ 0) with respect to the 
same set of histories C a . While not so difficult to satisfy as strong decoherence, it is 
still a very stiff requirement. For example, we shall see, in the linear oscillator models 
discussed below, that it is not very well satisfied there. However, if satisfied, medium 
strong decoherence would supply, in the case of an impure p, the same attractive features 
that medium decoherence yielded for the pure case, including the permanence of the past. 

Approximate medium decoherence for density matrices continues to be defined by the 
approximate satisfaction of the condition (foureight) and can be discussed for the type of 
model coarse grainings considered in this paper, where a fixed set of coordinates is dis- 
tinguished and the other coordinates included in the model are ignored. The mechanism 
of the formation of correlations between distinguished coordinates and ignored ones can 
continue to operate for each state in the density matrix. Now, however, there is the possi- 
bility for improvement in the effectiveness of approximate decoherence from summations 
over the states in the impure density matrix. 

To make this explicit, consider the Schrodinger evolution of an initial density matrix 
of the form (fourtwentythree b). The evolution described by eqs. (fourfifteen) - (fournine- 
teen) is as before except that each x acquires an index \i. In particular (fourtwenty) and 
(fourtwentyone) become 

D(a\a)= £ E (^'^E/^Cr X^ Wl ) , (fourtwentysix) 
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(fourtwenty seven) 
The additional sum over fj, can lead to further phase cancellations and more effective 
decoherence. We will see an illustration of this in the discussion of the oscillator models 
discussed below. 

(b) Decoherence in the Linear Oscillator Models 

Medium decoherence can be treated explicitly in the kind of oscillator models worked 
out by Feynman and Vernon [12] and Caldeira and Leggett [11], which we described in 
Section III. Generally these models assume an impure initial density matrix factored as in 
(eight), with the ignored coordinates representing a continuum of oscillators in a thermal 
state. 

An explicit illustration of the medium decoherence of a pure state can be found in 
the zero temperature limit of the linear oscillator model. In that limit all the ignored 
oscillators are in a pure state and the initial state of the distinguished oscillators may be 
taken to be pure. Denote the initial ground state of the bath by xo(Qo) an d the initial 
state of the distinguished oscillators by ^o(^o), so that the pure initial wave function of 
the whole system is ^o(xq)xq(Qq)- 

The Schrodinger evolution of this state can be represented in path integral form. At 
time t we have 

\I/ (x, Q, t) = l 5x5Q ex.p(iS[x(r) 7 Q(t)]/Ji)iPq(xq)xo(Qo) (fourtwenty eight) 

The integral is performed first over paths that start at xq and Qq at time £q and end at x 
and Q at time t, and then over xq and Qq. The integral over the distinguished coordinates 
x may be written as a sum over coarse grainings a\, ■ ■ ■ , a^ , % < t and an integral over 
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the paths restricted by the coarse graining, in the form 

»(*,«,«) = £ /fc««P(«t-(r),0(r)]/n)«,(x l) )x»(Oo). (/o^ent^) 

ai,— ,afc Ja 

This can be rewritten in the form 

/+oo 
dr S(x - r) x ak '" ai (r, Q) (fourthirty) 
-oo 

where 

X c " l '" ai (r,Q)= I SxSQ exp(iS[x(T),Q(T)]/h)i;o(xo)xo(Qo) (fourthirty one) 



ai, •••,<% " _:s 



and the integral is over paths, consistent with the coarse graining, that start at (xq, Qq) 
and end at (r,Q), including an integral over the values of xq and Qq. This is evidently 
the analog of (fournineteen) with r being a continuous index and 4>^(x) = S(x — r), branch 
independent and the same at each time. The overlap that occurs in (fourtwenty), of course, 
gives the decoherence functional. 

In the zero temperature limit, Caldeira and Leggett find for the imaginary part of the 
influence phase [cf. (6.8)] 

ImW\x , (T),x(T)]= - dt dt?£(t)ki(t-l/)£(tf) . (fourthirty two) 

4 Jq Jq 

Here, £(£) = x'(t) — x(t) and 

4M7 f n 

kj(r) = / duu cos(ut) , (fourthirty three) 

k Jq 

where O is the cutoff of the oscillator spectrum. This imaginary part of the influence 

phase favors contributions to the functional integral defining the decoherence functional 

(five) from values of £ near zero, or x'(t) close to x(t). 

A crude estimate of the time intervals by which coarse-grained alternatives of position 
{A^, fc } must be spaced in order to ensure decoherence may be obtained as follows: Approx- 
imate £(£) by a constant value d that is characteristic of the sizes of the intervals {A* }. 
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Find the time interval ^decoherence sucn that the integral in (fourthirtytwo) evaluated over 
that time interval is of order of magnitude unity. One then has a rough estimate of a time 
interval long enough for ImW/h to be large enough for exp(— XmW/h) to be small. The 
answer is 

^decoherence ~ ^ ex P f Jf^fi) ' (fourthirtyfour) 

assuming ^decoherence >> 1- Note that ^decoherence decreases as the coupling 7 is made 
stronger or the graining is made coarser or the number of oscillators becomes larger. In 
the present case of a pure initial state, the time scale for decoherence is essentially set by 
the cut-off, for fixed coupling and coarse graining. 

The improvement in decoherence from the sum over states in an impure density matrix 
is well illustrated in the linear oscillator models. In the high temperature Fokker-Planck 
limit, where many states contribute with nearly equal probability to the density matrix 
p^, the imaginary part of the influence phase is given by [cf. (ten)] 

rp 

lmW[x'(T),x(T)]= ^ — - / dt [x'it) -xit)] 2 . (fourthirtyfive) 

The large value of XmW suppresses contributions to the functional integral (five) defining 
the decoherence functional when x'(t) is significantly different from x(t), provided enough 
time has elapsed between successive alternatives so that significant values of XmW/h are 
built up, yielding approximate medium decoherence. Further, in this limit of very large 
Tg such that /cTg >> hfl, there is more efficient medium decoherence than is provided in 
the pure ground state example of (fourthirtyfour). 

Another feature that can be illustrated in the linear oscillator models is the permanence 
of the past. We have seen how the successive narrowing of the records implied by medium 
decoherence in the case of a pure state gives a natural explanation of the permanence of 
the past. For an impure state there is no satisfactory corresponding notion of generalized 
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record, but in the oscillator models it is still possible to show how the past becomes 
permanent. 

History is approximately permanent for a suitably restricted class of coarse grainings in 
the oscillator model in the Fokker-Planck approximation. The reason is that decoherence 
there is essentially local in time. More precisely, consider the following integral: 

/ Sx' / 5xM (x' k+v x k+ i) expli(S iYee [x'(T)] 
J a' J a "■ 

-5free[a;(r)] + W[x'(r), x(r)])/hJN (x' k , x k ) , {four thirty six) 

where the path integrals, as well as the integrals defining the actions and the influence 
phase, are over the interval of time from the time t k of a set of alternatives {cx k } to the 
time of the next set t k +\. The values x' k , x k are the endpoints at time t k and ^ +1 , x k+ \ 
are the endpoints at time i&+l- We consider any functions M and N. The expression 
(fourthirtysix) for the imaginary part of the influence phase implies that, if the intervals 
{A^, fe } by which the paths x(t) are coarse grained have a characteristic size d and the 
interval between t k and t k+ i is larger than the characteristic decoherence time scale [20a] 

i 2 

, (fourthirty seven) 



1 
^coherence ~ " 



n i 



v /2MkT B d 

then the "off-diagonal" terms in (foureight) will be very small and approximate decoherence 
will be achieved for the alternatives {ctfc} for a large class of functions M and N. 

Consider a set of coarse-grained histories defined by regions {A^, fc } at times t\, ■ ■ ■ , t n 
separated by time intervals longer than the characteristic decoherence time (fourthirty- 
seven). Fine-graining this set by adjoining further sets of intervals at similarly spaced 
times greater than t n does not affect the decoherence of those already present because 
the mechanism of decoherence exhibited by (fourthirtysix) is operative over a time scale 
of the order of ^decoherence about the time of each set of alternatives. Physically, that is 
reasonable. In this model, phases are carried away by interactions, local in time, of the 
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distinguished variables x a with the rest. Once dispersed among the continuum of oscil- 
lators described by the Q , they cannot be recovered by finer grainings beyond t n that 
involve the x a alone. To recover the phases, one would need a much finer graining that 
involved the whole set of variables. 

V. Distributions for Decoherence Punctionals 

As the review presented in the preceding two Sections makes clear, quantum theory can 
be organized into two parts: First, there is the calculus of amplitudes for histories or the 
bilinear combinations that are the decoherence functionals. The rules of this calculus derive 
ultimately from the principle of superposition. Second, there are the rules for deriving 
probabilities from these amplitudes, most generally, the several notions of decoherence of 
sets of alternative histories of a closed system. In this section, we show how the first part 
- the calculus of decoherence functionals for histories -- may be usefully reexpressed in 
terms of distribution functionals analogous to the Wigner distribution for alternatives at 
a single moment of time. 

In familiar quantum mechanics the probability that a determination of coordinates x a 
at one moment of time will yield a result in a volume V of the reduced configuration space 
spanned by these coordinates is 



PiY) = / dxp(x,x) , {tend) 

Jv 

where p(x', x) is the reduced density matrix on the reduced configuration space [cf. (four)] 
and dx is the reduced volume element. As is well known, the density matrix p(x , x) that 
gives such probabilities may be usefully expressed in terms of the Wigner distribution on 
phase space 

p {X + £/2, X - f /2) = / dPw{X, P)e^ P l n . (tenb) 
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The probability p(V) is then given by 



p(V) = dX J dPw (X, P) . {tenc) 



It is also true that the probability density for the momentum, P, conjugate to x is 
given by f dXw(X,P). In these properties w(X,P) is like a classical distribution on 
phase space. Other ways in which it is similar have been extensively discussed. (See e.g., 
[20, 21]) It differs from a classical distribution in that it is not in general positive and it 
does not provide analogs for all probabilities that are defined on classical phase space. For 
example, it would be incorrect to think of w(X, P) itself as a probability for a simultaneous 
determination of position and momentum. The calculus of amplitudes must therefore be 
supplemented by the rules that specify which classical quantities on phase space can be 
assigned probabilities. Ignoring those rules, and considering just the calculus of amplitudes 
for alternatives at a single time, we see from equations (tenb) and (tenc) how the calculus 
may be reformulated in terms of a distribution function on phase space, although not, in 
general, in terms of a positive one. In the following we shall give an analogous formulation 
for the calculus of amplitudes for time histories. 

The decoherence functional defined in equation (3.1) is the bilinear combination of 
amplitudes of which the diagonal elements give the probabilities of the individual histories 
in a decoherent set [cf (fourseven)]. In this sense, it plays the role for histories that the 
reduced density matrix does for alternatives at a single moment of time. We now construct 
a distribution functional for the decoherence function in much the same way in which the 
Wigner distribution was constructed from the reduced density matrix. 

We consider the partially coarse-grained decoherence functional D[x (t),x(t)] for his- 
tories that are fine-grained in the variables x distinguished by the initial coarse-graining 
discussed in the previous section. Introduce variables X(t) and £(t) that are the average 
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and difference respectively of the arguments of the decoherence functional: 

X(t) = -(x'(t) + x(t)) , (tenda) 

£(t) = x'(t) - x(t) . (tendb) 

In terms of these variables the fine-grained decoherence functional defined by equation 
(five) may be written 

D[X(t),£(t);X ,£o) = <*(£/) exp QU[X(t),£(t);X ,£o)) P (*o + f ,*0 - f ) , 

(tene) 
where 

A[X(t), £(t); X , ft) = S iTCC [X(r) + £(r)/2] - S ivcc [X(r) - £(r)/2] 

+W[X(t),Z(t);X ,Zo) . (tenf) 

A distribution G[R(t), X(t),Xq, ft) may be introduced for the decoherence functional 
by taking its functional Fourier transform with respect to £(£). We define G by the formula 

The functional G may be calculated from the decoherence functional by an inverse func- 
tional Fourier transform. Expressions like (teng) are to be interpreted as limits of mul- 
tiple integrals over paths that are piecewise linear between a discrete set of time slices 
to — 0, t\, • ■ • , tjy = T as the number of slices, N, tends to infinity. This is a standard way 
of defining path integrals; details for the particular integrals of interest are discussed more 
fully in Appendix (c). As an aid to the present discussion, however, it is useful to note 
that we always represent the integral such as that in the exponent of (teng) by discrete 
sums of the form 

rp 

J dT^(T)R(r) = ^ =1 e^ k R k , (tenh) 
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D[X(tU(t);X ,Zo) = JsRexp % - J dt^(t)R(t) 



G[R(t),X(t); X ,ft) • (teng) 



where e is the separation between time slices. The variable £o is thus not transformed and 
appears on both sides of the equation. 

The transform variable R(t) has the dimension of force. As we shall see later, 
G[X(r), R(t);Xq, £o) may be regarded as a classical distribution of R(t), given the path 
X(t) and Xq and £,q, although in general a non-positive one. It is therefore instructive to 
consider its moments. We define 

„ „,, , „,, w JtR R(h)-R(tn) G[R(T),X(T);X (h (, ) 

<m) ■••*<*•>>«= jmmWxW^M ' iten3) 

Clearly, < R(t\) ■ ■ ■ R{t n ) > c is a functional of the path X(r) and a function of Xq and £q- 
The first moment we define to be the average "total force" 

£{t, X , faX(T)] =< R(t) > c . (tenk) 

By "total force" we mean the force minus the inertial term, so that £ = is the effective 
or phenomenological classical equation of motion on the average. The deviations of the 
force R(t) from its expected value define the Langevin force C[t, Xq, £q; X(t)] for a given 
path X(t): 

C(t,X ,toX(T)] = R(t)-£(t,X ,toX(T)] . (ten/) 

The reason for these designations will become clear when we express the quantum-mechanical 
probability for an individual history a in a decoherent set of histories in terms of the dis- 
tribution G. 

The fully coarse-grained decoherence functional (3.5) is given by 

D(a',a) = [ Sx' [ SxD[X(t),C(t);X ,^) , (tenm) 

J a' J a 

where X and £ are connected to x' and x by (tend) and their range is restricted by 
the coarse-graining through (tend). The partially coarse-grained decoherence functional 
is given by (tene) and (tenf). If the further graining defined by the successive sets of 
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regions is coarse enough so that sufficient positive imaginary part of W is built up between 
one set of intervals and the next, then there will be a significant contribution to the 
integral defining the decoherence functional only for values of £(£) near zero and for a = 
a'. (See Figure 1). That is medium decoherence. Further, in the diagonal elements of 
the decoherence functional, which are the probabilities of the individual coarse-grained 
histories, the integral over £(£) may be carried out, to an excellent approximation, as 
though unrestricted by the coarse graining, provided the intervals are sufficiently coarse. 
(See Figure 1). When reexpressed in terms of the distribution G, the probabilities for the 
individual histories p(a) are, in this approximation, 

p(a) = J SX J SR J 6£ G[R(T),X(r);Xo,^)exp % - J dt^(t)R(t) 

where the integrals over £(£), including that over £q 5 are unconstrained. 

The expression (tenn) for the probabilities of the individual histories in a coarse-grained 
set has an especially transparent interpretation in the special case where the initial density 
matrix factors as in (eight) and a simple kinetic minus potential energy form (six) is 
assumed for the action Sf iee [x(r)]. A simplification following from factorization is that the 
influence phase W has no direct dependence on x'q and xq, and hence on Xq and £q> except 
through the paths x'(r) and x(t). Thus, we write simply W[x'(t), x(t)]. More accurately, 
in a time slicing implementation of the functional integral in which integrals over the paths 
are discretized as in (tenh), W is independent of Xq and £q. Further discussion is given in 
Appendix (b). When the action Sf Tee [x(r)] has a simple kinetic minus potential form the 
dependence of the rest of the exponent A[X(t),£(t);Xq,£q) may be easily isolated. An 
integration by parts in (tenf) yields 

A[X(t),£(t);X ,£o) = -£jMX + i[X(r),£(T)] , (tenna) 

where 

rp 

A[X(r)^(r)]= [ db\-tf(t)MX(t) - V (X(t) + £(i)/2) + V (X(t) - £00/2)' 
JO L 
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+W[X(t),£(t)]. (tennb) 

Then, under the ground rules for discretizing functionals discussed above, the functional 
A is independent of both £q and Xq. In particular, the only terms in the fine-grained 
decoherence functional (tene) that depend explicitly on £q are the density matrix and the 
surface term in (tenna). The integral over £q may therefore be carried out in (tenn), giving 

p{a) = f 5X f 5R f 5£ g[R(r),X(r)] w(X ,MX )exp % - f dt^(t)R(t) . (tennc) 

where w is the Wigner function defined by (tenb) and the distribution g[R(r),X(r)] is 
defined by 

<K£/)exp[^i[X(T),£(r)]] = JsRg[R(r),X(r)}exp[^J dt^(t)R(t)] . (tennd) 

These expressions acquire precise meaning in the time slicing implementation of the path 
integrals discussed in Appendix (b). The important point for the present discussion is 
that the assumptions that the initial p factorizes and that <Sf ree has a simple kinetic minus 
potential energy form lead to a factorization of the general distribution G into a distri- 
bution, w, of initial values Xq and MXq and a distribution, g, involving forces R(t) and 
paths X(t). Further, the equation of motion £(t, X(r)] and the Langevin force C(t, X(t)] 
become independent of the initial conditions and calculable just from the distribution g 
of force R(t) through expressions of the same form as (tenj), (tenk), and (tenl) with G 
replaced by g. These features allow a simple interpretation, which we shortly describe. 

The Wigner distribution by which the initial conditions are distributed in the expres- 
sion for the probabilities (tennc) is not generally positive and neither is the distribution 
g. However, the probabilities p(a) must be positive. It is not difficult to see how this 
comes about. Were the restriction on the range of integration arising from the coarse 
graining restored in expressions (tennb) or (tennc), the numbers p(a) they define would 
be manifestly positive. That is because they are expressions for a diagonal element of a 
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decoherence functional which is always positive [cf. (fourfive)]. This is not unlike other 
smearings of the Wigner distribution [44] which are known to give generally positive re- 
sults. Of course, the approximation in which the restrictions of the coarse graining on the 
integration are ignored may result in small negative probabilities, but, to the extent the 
approximation is good, these are equivalent to zero for physical purposes. 

The representation (tennc) allows the probabilities of decohering coarse-grained his- 
tories p(a) to be thought of as the probabilities of the histories of a system moving 
classically under the action of a stochastic force. The Wigner function gives the dis- 
tribution of initial conditions. The distribution g[R(r), X(t)] may be thought of as the 
distribution of total force R(t) acting on a system that describes the path X(t). Alterna- 
tively, if reexpressed in terms of £(t,X(r)] and £(t,X(r)] defined by (tenm) and (tenl), 
g £{t,X{r)} + £(t, X(t));X(t) may be thought of as the distribution of C(t) given the 
path X(t). In the approximation we have discussed, the unconstrained integration over 
£(£) in (tennc) leads to a functional 6- function that enforces the condition, R(t) = 0, that 
the "total force" on the system vanish. That is, it enforces the effective classical equation 
of motion, corrected by the Langevin force: 

£(t,X(r)]+C(t,X(r)] = . (tenne) 

In this equation the average total force £ (t, X(t)] is the known functional of X(r) defined 
by (tenk). The force £(£, X(t)] is distributed according to the distribution g. It is for this 
reason that we have called C the Langevin force; it can be thought of as noise. 

The probabilities p{a) for decohering coarse-grained histories are thus obtained from 
a mathematical description of classical dynamical system characterized by an equation of 
motion £(t,X(r)} together with distributed initial conditions and distributed noise. We 
should stress that this does not mean that quantum mechanics is equivalent to some kind 
of classical physics. For one thing, the distributions of noise and initial conditions are not 
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generally positive and that is certainly a non-classical feature. For another thing, there is 
the coarse graining needed for decoherence --an entirely non-classical requirement that 
must be satisfied before the p(a) may be considered as the probabilities of histories. 

The non-positivity that distinguishes quantum-mechanical distributions from classical 
ones may be regarded as the reason that Bell's classical inequalities [37] are violated in 
quantum mechanics, leading to important experimental tests of the theory. Of course, the 
Wigner distribution considered here is for continuous variables while Bell's discussion of the 
EPRB problem was for a discrete spin system. However, Feynman [32] showed that there is 
an analog of the Wigner distribution for spin-^ systems and that the departure of quantum 
mechanics from Bell's inequality is traceable to the fact that the analogous distribution is 
not generally positive. At the level of the calculus of amplitudes, the difference between 
classical and quantum mechanics is just the possibility of negative distributions. That 
possibility alone does not completely characterize the difference when we go beyond the 
calculus of amplitudes, because in quantum mechanics we also have the requirement of 
decoherence of histories. 

When the initial p does not factor, or when <Sf ree is not of simple kinetic minus potential 
energy form, the interpretation of the general expression (tenn) for the probabilities is 
less direct. However, building on the analogy of the special case, we may still think of 
G[R(t),X(t); Xq, £q) as a combined distribution of total force and initial conditions given 
the path X(r). Now, however, the distribution of the force R(t) is not independent of 
the initial conditions but depends on them. Further, the distribution of initial momenta 
is not given directly but only implicitly through the integral over £q- Finally, the equation 
of motion and noise depend on the initial values of Xq and £o- The integral over £(t) in 
(tenn) continues to enforce the classical condition that the total force R(t) = vanish. 
These features complicate the interpretation of (tenn) but they do not vitiate its validity. 

When the noise C(t, X(t)] is small compared the equation of motion term £ (t, X(t)] in 
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(tenne), we expect approximate classical determinism. More precisely we expect significant 
probabilities for histories correlated in time by the classical equation S(t,X(r)] = with 
small deviations produced by the noise. In the following Sections we shall analyse the 
circumstances where this is so. 

VI. Linear Systems 

(a) Equations of Motion 

We begin our discussion of the derivation of classical equations of motion for quantum 
systems by considering the simplest possible example -- linear systems. This is the class 
of models studied by Feynman and Vernon [12], Caldeira and Leggett [11], and many 
others and for which there is a wealth of information available on the specific forms of 
the influence functional, its dependence on the initial condition, etc. Either implicitly or 
explicitly, equations of motion have been considered for these models by several authors. 
We are thus on familiar territory. 

Linear systems may be characterized precisely, following Feynman and Vernon, by the 
following two requirements: 

1. A free action for the variables distinguished by the coarse graining with a kinetic energy 
that is quadratic in the coordinates and their velocities 

rT pi i , 

S bee [x(t)] = dt - x\t)Mx(t) - - x\t)Kx(t) . (eleven) 

J L2 2 J 

We use here and throughout an obvious vector notation so that M and K are positive 

constant matrices and (l/2)(x<Mx) = (l/2)T, a i ) x a M a i ) x p , etc., in this case where the 

variables are real. 

2. An influence phase that is at most quadratic in the variables x(t). Its most general 
form has been deduced by Feynman and Vernon [12] from general symmetries and 
quantum-mechanical causality. (See also the exposition in [22].) It consists of terms 
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linear in the distinguished variables x(t) and x (t) and terms quadratic in them. The 
linear terms may be eliminated by a time-dependent shift in x(t) and for simplicity we 
imagine this has been done. The general form is then 

W \x'(t),x(t)] =\ I dt [ dt' \x'(t) -x{t)} ] \k(t,t')x'(t')+k*(t,t')x(t')] 
2 Jo Jo 

= - I dt [ dt'[x'(t) -x(t)]Uk R (t,t')[x'(t') + x(t')] + ik I (t,t')[x'(t') -x(t')]\ , 

(twelve) 
where k(t,t) is a complex matrix kernel with real and imaginary parts k R (t,t) and 
kj(t, t') respectively. We shall explicitly assume that W depends on xq and x'q only 
implicitly through its dependence on the paths x'(t) and x(t) and not explicitly as in 
the general case. 

Only under very restrictive conditions will the influence phase be exactly quadratic 
as in (twelve) with no explicit dependence on x'q and xq. Such an influence phase will 
certainly follow if (i) the action of the ignored variables, So[Q(t)], is quadratic in the 
Q's; (ii) the interaction between the distinguished and ignored variables is exactly linear 
in each, giving 

rT 

S- mt [x(t),Q(t)} = - dtx^(t)f(Q(t)) , (thirteen) 

Jo 

where the force / is homogeneous and linear in the Q's; and (iii) the initial density matrix 
factors as in (eight) with a density matrix pb(Q'qi Qo) that is of the form 

Pb(Q'oiQo) =exp [-B(Q' ,Qo)] , (thirteena) 

where B(Q'q,Qq) is quadratic in its arguments. Under these conditions the integral in 
(three) is a Gaussian and a quadratic influence phase results. The models defined by 
(eleven), (twelve), (six), and (thirteena) are, in fact, just those considered by Feynman 
and Vernon [12] and Caldeira and Leggett [11]. 
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A useful explicit example, studied in [11, 12], is the case of a single coordinate, x, 
interacting with an assembly of "bath" oscillators, with an initial condition that factors as 
in (eight) and with a state of thermal equilibrium for the "bath" oscillators at temperature 
T = (k(3) . The oscillators are described by a free action 



So[Q(t)] = J2j dt \ m {Q A ^) 2 -\mu\(Q A (t) 
The force / entering (thirteen) may be written explicitly as 



(thirteenb) 



A 



A 



(thirteenc) 



with coupling constants C A , and the function B is independent of x'q and xq and given, 
up to an additive constant (determinable through normalizing p_g), by the expression 



B {Q'o, Qo) = E 



mcc^4 



A 



2h smh(h(3u>A) 



(Qo 4 ) 2 + (Qo) 2 ] co S h(h(3u; A ) - 2Q' A Q$ 

(thirteend) 



Then, as in [11] or [12], 



C 2 
k R (t,t') = -J2 —^ smu A (t-t') 
A muJ A 



and 



C 2 /l \ 

kj(t, t') = Y^ — ^- coth I -hj3u A cosu A (t - t') 
^ mu A \2 J 



(thirteenea) 



(thirteeneb) 



The influence phase quoted in eq. (ten) for a distinguished oscillator interacting with a 
high temperature thermal bath is a special case of these expressions in which a continuum 
limit of oscillators with the special couplings 



p D (u)C 2 (u) =4mM-yu 2 /7v 



(thirteen/) 



below the cutoff was assumed, where pjj is the density of states. With this coupling, we 
have the following, with the high temperature limit needed only in (6.10b) and the high 
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cut-off, which facilitates phase dispersal, used in both equations: 

k R (t, t') = -4M~/5'(t - t') , (fourteena) 

kj(t 7 t') = 6(t — t) , (fourteenb) 

n 

where the S function occurring in retarded integrals is a distribution such that 

/ g(t') [S'(t - t')] dt' = \ g'(t) . (fourteenc) 

J— 00 ^ 

[Eq. (fourteen a) gives the renormalized value of k R (t,t f ). Infinite terms in (twelve) pro- 
portional to 5(t — t) that arise from the continuum limit taken by Caldeira and Leggett 
have been absorbed in a renormalization of the frequency of the distinguished oscillator in 
its free action.] Although the influence phase is exactly quadratic only under restrictive 
circumstances, the linear cases supply useful models for more general ones, as we shall see. 

With these preliminaries, we can now give a derivation of the classical equation of 
motion for these linear systems. The equation we shall derive is, of course, the same as 
considered, for example, by Caldeira and Leggett. Even in this linear case, however, we 
believe that there are several important new features of this derivation. It is consistent 
with the general discussion of the average equation of motion and noise in Section V. The 
probabilities that coarse-grained histories of the distinguished particle are correlated in 
time by equations of motion are explicitly considered. The form of the phenomenological 
equation of motion is derived from a consideration of these probabilities. The amount of 
coarse graining necessary for the decoherence of histories and their classical behavior is 
discussed quantitatively, and the connection between decoherence and quantum noise is 
made explicit. Most importantly, the derivation suggests how the generalization to the 
nonlinear case is to be carried out. 

The imaginary term in the influence phase (twelve) gives rise to decoherence between 
the trajectories of the distinguished variables, provided that the coarse graining is such that 
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the integral of this term is sufficiently large for different coarse-grained histories so that 
the corresponding "off-diagonal" elements of the decoherence functional are exponentially 
small. To exhibit this decoherence explicitly it is useful to change variables, in the integral 
(five) defining the decoherence functional, from x (t) and x(t) to the average and difference 
X(t) and £(t) defined by (tend). The exponent in the decoherence functional (five) can be 
explicitly expressed in terms of X(t) and £(£) using (twelve). Denoting this exponent by 
A[X(t),£(t)] as in (tenf), one finds after a few integrations by parts that 

D(a',a)= ! 6X6£6(Zf) exp(iA[X(rU(r)]/n)p(X + ^/2,X -^/2) , 

J (a', a) 

(fifteen) 
where 

A[X(t),£(t)] = -^MX Q + f dt^(t)e(t) + - f dt f dt^Ht)kj(t,t f )C(t') . 

Jo 4 JO JO 

(sixteena) 
Here, e(t) is the average "total force" 

e(t) = -MX(t) - KX(t) + [ dt'k R {t,t')X(t') , (sixteenb) 

Jo 

and kj is defined to be symmetric in its argument, so that the limit of integration on t 
can be extended to T. The quantities Xq, Xq, etc. are understood in the usual path 
integral sense as finite difference expressions in a time-sliced implementation of the path 
integral. (See Appendix (b).) The integrals over X(t) and £(£) are constrained by the 
coarse graining defining the histories a' and a. 

The imaginary term in (6.12a) leads to decoherence. We shall give below, in the 
discussion of noise, an explicit construction of the kernel kj that shows it to be a positive 
kernel. The imaginary part of W is proportional to kj and occurs in the expression 
ex.p(iW/h), giving a decreasing exponential. If the graining defined by the successive sets 
of regions {A^ fe } is coarse enough so that sufficient positive imaginary part of W is built 
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up between one set of intervals and the next, then there will be a significant contribution 
to the integral defining the decoherence functional only for values of £(t) near zero and for 
a = a'. (See Figure 1). That is medium decoherence. 

If only values of £(£) near zero contribute significantly to the integral (fifteen), then 
in the diagonal elements of the decoherence functional, which are the probabilities of 
the individual coarse-grained histories, the integral over £(£) may be carried out, to an 
excellent approximation, as though unrestricted by the coarse graining when the intervals 
are sufficiently coarse. (See Figure 1). The integral over £q leads to the Wigner distribution 
(tenb) as in (tennc). The result of the unrestricted Gaussian integrals over the rest of £(t) 
is again a Gaussian functional. We thus obtain for the probabilities of the individual 
histories in the coarse-grained set, the following expression : 

p(a) ^ / SX [de*(fc//47r)l 5 exp[-i/" dt f dt'e\t)kf v (t, t')e(t')] w(X , MX ) , 
J a L J L n Jo Jo J 

(eighteen) 
where k l j lv is the inverse kernel to kj. The integral is over all paths that proceed from 
t = to t = T and lie in the class corresponding to the coarse-grained history a. The 
integral includes an integration over the initial and final endpoints Xq and X j respectively. 
Again, for further details and an explicit representation of (eighteen), see Appendix (b). 

The Gaussian exponential in (eighteen) means that, for given Xq and MXq, the histo- 
ries with the largest probabilities will be those with e(t) = 0, that is, those for which the 
time evolution is correlated according to the effective average classical equation of motion 

e(t) = -MX(t) - KXit) + / dt'k R (t,t')X(t') = . (nineteen) 

Jo 

This is, of course, not the equation of motion following from the free action of the dis- 
tinguished oscillator. It differs by the additional force, non-local in time, that arises from 



Dowker and Halliwell [40] have obtained analogous expressions in linear models for the 
probabilities of histories defined by a finite number of "Gaussian slits" . 
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the interactions of the distinguished variables with the rest. The presence of such a force 
will in general mean that energy is not conserved, leading sometimes to dissipation. Al- 
though non-local in time, the additional force in (nineteen) is retarded, expressing classical 
causality. The origin of this retardation can be traced to the retarded form of the general 
influence phase (twelve). That, in turn, follows from quantum-mechanical causality the 
fact that the decoherence functional has a trace in the future and a density matrix differ- 
ing from the unit matrix in the past. Causality in quantum mechanics thus implies the 
causality of classical physics. 

A special case of a linear system is the Fokker-Planck limit of the oscillator model, for 
which the influence phase is exhibited in (ten). With the corresponding kji of eq. (6.10a), 
the equation of motion away from t = becomes 

e{t) = -MX(t) - KX{t) - 2M 1 X{t) = . (twenty) 

This is local in time, but that is a special property of the way in which the limit of a 
continuous spectrum of oscillators was taken in the Caldeira-Leggett model, not a general 
one. In that limit, eq. (twenty) explicitly exhibits the familiar form of frictional dissipation, 
not necessarily a general characteristic of the additional force in (nineteen) . 

The individual classical histories in (eighteen) are distributed according to the probabil- 
ities of their initial conditions Xq and Pq = MXq given by the Wigner function w(Xq, Pq). 
Although the Wigner function is not generally positive, we know, as discussed in Section 
V, that apart perhaps from small errors introduced by the approximation in which the 
constraints of the coarse graining on the ^-integrations were neglected, the result of the 
integral (eighteen) must be positive even though the Wigner function is not. 



We are thus for simplicity ignoring the terms proportional to x(0) that arise when the 
integral in (sixteen b) is carried out using (fourteen). For further discussion see [36]. We 
thank J. P. Paz for a discussion of this point. 

For explicit examples of this see [45] . 
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(b) Noise and Predictability 

The distribution of probabilities for histories (eighteen) predicts the largest probability 
for histories obeying the classical equations of motion but also predicts probabilities for 
deviations from classical predictability. Those give the noise, including quantum noise. The 
same interactions (of the distinguished variables with the others) that carry away phase 
information to produce decoherence also produce the quantum and classical-statistical 
buffeting of the trajectory of the distinguished variables that constitutes the noise. 

In Section V we showed how the probabilities of coarse-grained histories p(a) could 
be thought of as the probabilities of a classical dynamical problem with (generally non- 
positive) distributions of force and initial conditions. We now specialize that discussion to 
the linear systems of the present Section. 

Since we have assumed factorization of the initial p, we may compute separate distri- 
butions for the initial conditions and for the total force as in (tennc) . The initial conditions 
are distributed according to the Wigner distribution as shown by that equation or directly 
from (thirteen). The total force is then distributed according to the distribution function 
g[R(r),X(r)] defined by (tennd). Using (6.12a) for A and (tenna) to define A, we may 
calculate g[R(r),X(r)] directly. It is 

g[R(r),X(r)} = [det(k I /2iT)}- l 2 



x exp 



rp rp 

I f dt [ dt'(R(t) - e(t))hf Y (t,t')(R(t') - e(t')) 

hJo Jo 



(twenty a) 



where e(t) is given by (6.12b) and the precise meaning of the inverse kernel k l j nv (t,t') is 
discussed in the Appendix. 

Evidently, we have, for this case, 

£(t, X(t)} = (R(t)) c = e(t) . (twentyb) 
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The equation of motion defined as the expected value of R in Section Vb, therefore, co- 
incides with the equation of motion e(t) whose correlations are favored by the probabili- 
ties (eighteen). The Langevin force C(t,X(r)] that governs the deviations from classical 
predicability is distributed according to g[e(t) + C(t,X(r)],X(r)]. As is easily seen from 
(twentya) for the linear models under consideration this noise is distributed with a positive 
Gaussian probability distribution that is independent of the path X(r). To emphasize this 
we write 

C(t,X(r)]=£(t) 



in the linear case, and £(t) is then distributed according to 



1 



T 



T 



, , dt I dt't\t)kf Y (t,t'y(t') 

ft./o Jo 



[det(fcj/27r)] 2 exp 
The spectrum of the Gaussian noise is summarized by the formula 



(twentyc) 



(twenty d) 



;*(t)*(0) c = (ft/2)*j(t,t')- 



(twenty five) 



In the Fokker-Planck limit of the oscillator model, we have 



'^(tW^^AM^kTB^t-t 1 ) 



(twentysix) 



giving rise to the model for Brownian motion discussed by Caldeira and Leggett[ll]. 

For linear systems instructive, explicit expressions for the functions kji(t,t'), and 
kj(t,t'), which describe the influence of the ignored variables on the distinguished ones, 
may be obtained in terms of quantum- mechanical expectation values of the force f(Q) 
defined by (thirteen) . As these are straightforwardly derived as special cases of the similar 
formulae, applicable to non-linear situations, to be discussed in Section VI, we shall just 
quote the results here. 

Consider the Hilbert space of square-integrable functions in the ignored Q's. Define 
an expected value, ( )q, of an operator A(t) evolving by the Hamiltonian of Sq by 



(A(t)) = Sp[A(t) PB } 
53 



(twentyseven) 



where Sp denotes the trace over the ignored variables. The subscript zero means that 
the time dependence of the operators inside the expected value is calculated using the 
Hamiltonian Hq of the Q's alone neglecting interactions with the x's. 

As a consequence of our convention that the influence phase has no terms linear in 
x(t), the expected value of f(Q(t)) vanishes: 



f(Q(t))) =0. (twentynine) 

The real and imaginary parts of the kernel k(t, t) may be expressed in terms of expected 
values of fluctuations in the force as follows 

hk R (t,t') = i([f(Q(t)), f(Q(t'))]) Q , {thirty) 

hk^t') = ({f(Q(t))J(Q(t'))}) Q , {thirtyone) 

where [ , ] and { , } denote the commutator and anticommutator respectively and the 
matrix elements of k are understood to be the tensor product of the /'s. Expression 
(thirtyone) shows explicitly that the kernel kj(t,t) is positive in the sense that 



r r 

dt dt'^{t)kj{t,t')^{t')>0 (thirtytwo) 

JO JO 



for any real vector £(£). The same equation and (twentyfive) demonstrate that the spectrum 
of the random Gaussian force in the Langevin equation e{t) + £(t) = is directly given by 
the quantum correlation function of the fluctuation in the force f(Q), viz. 

(£(t)£(t')) c = l({f(Q(t))J(Q(t'))}) Q . (thzrtythree) 

Expressions (thirty) and (thirtyone) lead to the essential content of the fluctuation- 
dissipation theorem. If pg is diagonal in the energy representation defined by the Hamilto- 
nian of Sq as it is for equilibrium distributions, then kji(t, t') and kj(t, t') are functions of 
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the time difference t — t , and their spectral weights are simply related. More specifically, 

it follows from the symmetries of commutator and anticommutator that we could write 

/•oo 

kR(t,t) = / dw kji(u) sin uj(t — t') , (thirty four a) 

JO 

/•OO 

kj(t,t) = / (ko kj(oj) cosuj(t — t ) . (thirty f ourb) 

JO 

li Pi is the probability of an energy eigenstate \i) with eigenvalue E( in the density matrix 
PB, we have 

hk R (u) = ^2 i(Pi ~Pj) /(Q(°)) 3/ $( E j ~ E i ~ huJ ) ' (thirty fivea) 
ij 

hk^u) = J2 (Pi+Pj)\(i\f(Q(0))\j)\ S(Ej -Ei- hu) . (thirtyfiveb) 

ij 

Then, in a thermal bath where pi = exp(—[3Ei)/Z, we recover the famous relation (see, 
e.g., [14]) 

hkj(uj) = h coth I — — J kji(u) . (thirtysix) 

This connects the kernel k R (t, t') governing the effective force in (nineteen) with the kernel 
kl(t,t) governing the fluctuations in (twentyfive) . This connection is the fluctuation- 
dissipation theorem. 

An important fact that emerges clearly from these linear models is that the same 
coarse graining and interactions that accomplish decoherence also lead to dissipation and 
noise. The fluctuation-dissipation theorem derived above is a well known example of the 
connection between noise and dissipation. The connection between decoherence and these 
two phenomena appears to have been less widely stressed. We now consider it. 

There is, in effect, a competition between decoherence and classical predictability. Con- 
sider, for example, a model of an oscillator interacting with a thermal bath. Increase the 
temperature of the thermal bath, Tg, and decoherence is more effective. The characteristic 
time td eco h, by which successive intervals of typical length d must be spaced to give de- 
coherence, decreases with Tg according to (fourseven). However, deviations from classical 
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predictability expressed by (eighteen) also increase. To ensure both classical predictability 
and decoherence we must consider a further limit, the limit of high inertia. In the present 
model, that is the limit of large M, and the exponent in (eighteen) can be written 



M f T ,. N-v ,„ lT „ -12 



dt 



X + M~ L V'{x) + 27X . (thirtyseven) 



8-fkT B J 

Assuming that V does not increase faster than M, we see that the probabilities for his- 
tories will become sharply peaked about the certainties implied by a classical equation of 
motion in the limit of large M/Tg even as Tg itself is becoming large to ensure efficient 
decoherence. 

Of course, in realistic situations the parameters of a given system, such as the mass, 
are fixed. The same kind of limit of high inertia can be achieved, however, by considering 
coarser and coarser graining of a kind that keeps increasing the inertia of the variables 
distinguished by the coarse graining. Take the case of coarse grainings defined by hy- 
drodynamic variables that are integrals, over suitable volumes, of densities of exactly or 
approximately conserved quantities such as mass, energy, or momentum. By making the 
size of these volumes larger, the resistance to noise can be increased. In the present model, 
decoherence and classical predictibility can be achieved only by varying the parameters of 
the model. In realistic situations they are achieved by a suitably coarse graining. 

Classical behavior requires sufficient coarse graining and interactions for decoherence 
but sufficient "inertia" to resist the deviation from predictability that the coarse grain- 
ing and interactions produce. Traditionally other descriptions have been given of the 
requirements for classical behavior of measured subsystems. Large action or high quan- 
tum numbers are often mentioned. While such criteria are not as precise or as complete 
as those deduced here, it can be seen from simple dimensional arguments that in typical 
situations an action will be large compared to % when the two requirements of decoherence 
and sufficient inertia are satisfied. Let us consider a one-dimensional oscillator model of the 

56 



kind just discussed in the high temperature limit. Let £ dyn be the shortest dynamical time 
scale of interest and assume that the coarse graining is characterized by sets of intervals 
of characteristic size A separated in time by i d yn- From (ten), it follows that decoherence 

requires 

M~fkT B 



ft 



(t dyn A 2 ) » 1 . 



(sixthirtytwo) 



From (thirtyseven), the requirement of sufficient inertia is 

2 



M 



*dv 



7 JfeT B ~ dyn \ t 2 



» 1 



(sixthirty three) 



dyn, 



if we assume t dyn « t rc i axat ion = 1 /l- 

These two requirements may be reexpressed in terms of the characteristic scale of 
classical actions 



S ~ t dyn M 



t d 



(sixthirtyfour) 



yn 



and the thermal correlation time ^thermal — fc/kTg. One finds 



(S/h) » (I/77) 



(s/h) » v 



(sixthirtyfive) 



for the requirements (sixthirtytwo) and (sixthirtythree) respectively, where rj is the ratio 



V ^dyn/ ''thermal J V' : dyn/' : 



yn/ ^thermalj ^dyn/ ''relaxation 



)■ 



(sixthirty six) 



Whatever the size of rj, the relations (sixthirtyfive) imply that S/h >> 1 in the classical 
limit. 

VII. Equations of Motion for Non-Linear Systems 



In this Section we generalize the results of the preceding two Sections to the more 
realistic situation where the action Sf ree [x(r)], and the influence phase W[x'{r), x{r)\ are 
not necessarily quadratic functionals of their arguments. We begin in (a) by deriving some 
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useful general properties of the influence phase. These are used in (b) to derive the form 
of the phenomenological equations of motion and analyze the restrictions on the coarse 
graining that permit the histories to stay close to solutions of those equations with high 
probability. In (c) we derive the classical causality of these equations of motion from 
quantum-mechanical causality. The linear theory, discussed in the preceding two Sections, 
is recovered in (d). 

(a) General Relations for the Influence Phase 

The influence phase W[x'(t),x(t)] is defined by the functional integral (3.5) over the 
ignored variables Q . We continue to assume a factored initial condition as in (eight), 
so that W has no explicit dependence on Xq and xq. A useful operator expression for 
W[x (r), x(t)] may be derived by noting the following: The integrals over the Q(t) in (3.3) 
are over paths between t — and t = T that are unrestricted except at their initial and 
final endpoints. They may, therefore, be thought of as defining the unitary evolution of a 
family of operators p(x'q, xq) in the Hilbert space TCq of square-integrable functions of the 
Q's. The dynamics of this evolution of the Q's is specified by the action 

S q [x(t),Q(t)} = S [Q(r)]+S int [x(r),Q(r)} , (thirtyeight) 

which depends on the path x(t). There is a corresponding Hamiltonian operator on 'Kn- 
it we assume that the interaction is local in time, specifically such that 

-T 
'0 
then that Hamiltonian at time t depends only on the instantaneous value of x(t), viz.: 



Sint[ x ( T )iQ( T )]= dt L- mt (x(t),Q(t)) , (thirtynine) 

JO 



H Q (x(t)) = H Q + H int (x(t)) . (forty) 

Here Hq is the Hamiltonian of the Q variables, omitting their interaction with the x 
variables, corresponding to the action 5*0 [Q]. The operator effecting the unitary evolution 
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generated by this Hamiltonian between times t and t is 

r * f i 

u t",A x ( T )] =T ex P[-^ / dt H Q (x(t))\ , (fortyone) 

where T denotes the time ordered product. 

Write pq for the density operation on Tin whose matrix elements are 

(Q'o \pb\ Qo) = Pb(Q'o, Q6) > (fortytwo) 

where Pb(Q'o^Qo) is the factor of the initial density matrix (eight) referring to the Q's. 
Utilizing the C/'s defined by (fortyone) and the pg defined by (fortytwo) the path integral 
relation (three) defining the influence phase W[x'(r)] may be reexpressed [see Appendix 
(d)] as 

exp(z W [x'(t),x(t)] /n) = S P ^U Tfi [x'(T)]p B ul {x(T)]} , (fortythree) 

where Sp denotes the trace operation on the Hilbert space Hq. We shall now use this 
relation to derive some useful general properties of W[x (r), x(t);xq, xq). 

First, it is an immediate consequence of Sp(A*) = (Sp(A))* and the form of right-hand 
side of (fortythree) that interchanging x'(t) and x(t) on the left-hand side is equivalent to 
complex conjugation. (This elementary result also follows directly from (three), as noted 
in [12, 22].) Thus, 

IZe W [x (r),x(r)] = — IZe W [x(t),x (t)] , (forty four a) 

Xm W [x (t),x(t)] = -\-Xm W [x(r),x (r)] . (fortyfourb) 

In particular, if the influence phase is written as a functional W [X(r), £(t)] of the average 
of x'(t) and x(t) and the difference between them [eq. (tend)], then lZe(W) is an odd 
functional of £(t) while Xm(W) is an even functional. 
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As shown by Brun [27] in the following paper, an elementary application of Schwarz's 
inequality shows that 

exp(-JmW[X(r),e(r)]) < J2 t P^^ U Tfi^ T ^ U T,oW(r)]\A)\ < £.W = 1 > 

(forty four a) 
where p^ are the eigenvalues and \ipj) the eigenvectors of pg. Thus XmW[X(r)j £(r)] 
is positive, which is essential for the convergence of the functional integral defining the 
decoherence functional as well as decoherence itself. 

The expression (fortythree) may be used to find convenient operator expressions for 
the coefficients of the expansion of W[X(t), £(t)] in powers of £(£). Generally, 



W[x{t),Z(t)] = W[x{t),0]+ J dtt\t)(jj^y 

+ \ f dt [ dt't\t)( JYt(*r) ^') + ---- (forty five) 

2 Jo JO \W)<W)/f(t)=0 

For the leading term in (fortyfive), we have, evaluating the right-hand side of (fortythree) 
at £(t) = x'(t) - x(t) = 0, 

exp(W [X(t), 0] /%) = Sp{u Tfi [X(t)} Pb U^ q [X(t)}) . (fortysixa) 

Using the cyclic property of Sp and the unitarity of Uj> q[X{t)\ it is easy to see that the 
right-hand side of (fortysix) is unity. Thus, the leading term in (fortyfive) vanishes: 

W[X(r),0] =0 . (forty sixb) 

To evaluate the next term in the expansion (fortyfive), we must consider the derivatives 
SU Tfi [X(t) ±f(i)/2]/<5f(i). To do this we introduce the definition 

The operator F(x(t)) is an operator in the Schrodinger picture corresponding to the unitary 
evolution operators we have been using. It is a function of x because L- m ^ is a function of 
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x and it becomes a function of t because x is a function of t. It is an operator representing 
the force on the distinguished coordinates x(t) due to their interaction with the rest of the 
system. 

Carrying out the indicated differentiations of the t/'s yields 

(sU Tfi [X(T)±Z(T)/2]/5at)) m=0 = ±(i/2WT,t [X(t)] F(X(t)) U tfi [X(r)] 

= ±(i/2h)F(t,X(r)}. (fifty) 

The operator F(t, X(t)] is the Heisenberg picture representative of the Schrodinger picture 
operator (fortynine). It is a function oft but also a functional of the path X(r). We indicate 
this dependence by writing f(t, X(r)], using a parenthesis on the left to indicate that it is 
a function and the bracket on the right to indicate that it is also a functional. 

With the result (fifty) , it is only a short calculation to find the coefficient of the linear 
term in (fortyfive). It is 

(SW/6£(t)) m =0 = S P{F(t, X(r)\ Pb \ . (fiftyone) 

If we define the Heisenberg picture expected value by 

(A) = Sp(A PB ) , (fiftytwo) 

eq. (fiftyone) may be written in the compact form 

(8W/6Z(t)) m=0 = (F(t, X(r)]) . (fiftythree) 

As required by (fortyfour), this contribution to the part of W that is odd in £(t) is purely 
real. 

The coefficient (5 W/5£(t)S£(t'))fn\ = Q of the quadratic term in the expansion of W in 
powers of £(t) is similarly evaluated. The expansion of the Ws on the right-hand side of 
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(fortythree) will result in products of first derivatives, such as those in (fifty), but also in 
second derivatives. Those second derivatives yield a term in the expansion proportional to 



S(t - t') Sp ( [< t (t, X(t)}, p B ] ) , (fifty four 



where H^At, X(t)\ is the Heisenberg picture representative of the operator d H- m ^/dx(t) . 
A term like (fiftyfour) vanishes, of course, for a linear interaction like (thirteen) . It formally 
vanishes in the non-linear case too because it is the trace of a commutator. Of course, 
that is a delicate issue in the case of unbounded operators, as the non-vanishing value 
of Sp( [IT, Q] ) shows. However, in the present case, where p^ is bounded and H[ nt is a 
function of the Q's and not of their conjugate momenta, we may reasonably assume that 
(fiftyfour) vanishes and we shall do so in what follows. 

The remaining contribution to the coefficient of the quadratic term in the expansion of 
W[X(t),£(t)] comes from products of first derivatives of t/'s such as those in (fifty) and 
is straightforwardly evaluated. One finds 

(5 2 W/5am(t%(t)=0 = (i/2ft)({AF(i, X(t)], AF(t',X(r)}}) , (fifty five) 

where {, } denotes the anticommutator and AF is the operator 

AF(t, X(t)] = F(t, X(t)] - (F(t, X(t)]) (fiftysix) 

representing fluctuations in the force F about its mean. We note that, as required by 
(fortyfour), this contribution to the even part of W[X(r), £(r)] in £(t) is purely imaginary. 
When divided by i, it is also manifestly positive in the sense of (thirtytwo). 

With these preliminaries we may now derive the non-linear equations of motion and 
discuss their form. 

(b) Non-Linear Equations of Motion 
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We consider a set of alternative coarse-grained histories specified at a sequence at 
times ti, • ■ ■ , t n by sets of exhaustive and exclusive regions of the x's which we denote by 
{A* }, {A| 2 }, ■ • ■, {A" n }. The decoherence functional for such sets of histories is given 
by (five). We assume that the regions and times are chosen so that there is a negligible 
contribution to the path integrals in the decoherence functional except when £(£) = x (t) — 
x(t) is small. We expect to have such coarse grainings, for example, if the imaginary part 
of the influence phase W[X(r), £(r)] has its minimum at £(£) = 0. Our assumption about 
the integrals in (five) implies the decoherence of a set of alternative coarse-grained histories 
(cf. Figure 1) and the following formula for their probabilities: 

p(a)<* [ 6X5Z8(Zf) exp{z(S frcc [X(r) + £(T)/2] 

- S iiec [X(t) - £(r)/2] + W [X(r)^(T)])/n)p (X + Co/2, X - ^i^ty seven) 



The functional integral is over paths in both X and £ as restricted by the coarse-grained 
history a. 

If only small values of £(£) contribute to the integrals in (fiftyseven) , we may make a 
further approximation by expanding the exponent in powers of £(£), utilizing the expansion 
in (fortyfive), up to the quadratic terms. We then have an integral for the probability of 
a history that is of precisely the same form as the one occurring in the discussion of the 
linear theory, (eighteen), with e(t) replaced by 

S{t) = 6Si Txft) T)] + (F(t ' X{T)]) (fiftyeight) 

and with kj(t,t f ) replaced by 

tfj(M') =h- 1 ({AF(t,X(r)lAF(t',X(T)]}) . (fiftynine) 



In writing out these identifications we have made use of (fortysix b), (fiftythree), and 
(fiftyfive) for the expansion coefficients of W[-X"(t),£(t)]. Eq. (fiftynine) shows the kernel, 
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Kj(t,t ;X(t)] is manifestly positive -- a necessary condition for the mechanism of deco- 
herence being discussed, not to mention the convergence of the integral (fiftyseven) with 
the expanded exponent. 

Under the assumption that only a narrow range of £(£) near zero contributes to the 
integral (fiftyseven), it is a good further approximation to neglect the constraints on the 
ingegration range of £(£) arising from the coarse graining. The resulting Gaussian inte- 
grals can then be carried out, yielding an expression for p(a) that is the generalization of 
(eighteen) : 

p(a)= / 5X [det(Kj/27r)]"5 

Ja 

rp rp 

- [ dt [ dtf S ] (t) Kf v (t,t f ;X(T)]£(t')} w(X ,P ) , (sixty) 
n JO JO J 



x exp 

L n 

where -ftTV nv (t, t ,X(t)] is the inverse kernel to Kj(t,t ;X(r)] and Pq is the momentum of 

the free action expressed in terms of Xq and Xq. The measure 5X is discussed in the 

Appendix (b). 

The derivation and analysis of the equations of motion now proceed as in the linear 
example, with important differences that we shall mention. The Gaussian form of the 
exponent in (sixty) means that for given Xq and Pq the most significant contribution 
comes from the histories with £{t) ~ 0, that is, those whose evolution in time nearly 
follows the effective classical equation of motion 

£(t) = 6SiT ^f t) (T)] + (F(t, X(r)]) = . (sixtyone) 

The probabilities predicted by (sixty) are, therefore, those of an ensemble of classical 
histories individually correlated in time by the equation of motion (sixtyone) and with 
initial conditions distributed according to the Wigner function w(Xq, Pq). 

The first term in (sixtyone) is the equation of motion of the distinguished coordinates 
x a in the absence of any interaction with the remaining coordinates Q . The second term 
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is expected value of the force arising from that interaction. This is a functional of the tra- 
jectory and will, in general, be non-local in time. As we shall show below, it is retarded as 
a consequence of quantum-mechanical causality. It typically leads to dissipation, although 
under some conditions the energy might actually increase. The phenomenological force 
is also generally dependent on the initial condition pg through (fiftyone). The familiar, 
phenomenological equations describing, say, dissipative friction, are characterized by a few 
parameters independent of initial conditions and are the result of further approximations to 
(sixtyone). These are typically good in situations where there is a significant contribution 
only from retarded times that are short compared to the relaxation times of that part of 
the bath that interacts significantly. If the distinguished system has energy large compared 
to /cTg, it will lose energy to the bath on the average. The result is a dissipative phe- 
nomenological equation, local in time, with parameters independent of initial conditions, 
like (twenty). 

Viewed as a generalization of the linear case, the important point about the equation of 
motion (sixtyone) is that both the free part of the equation of motion and the contribution 
from the interaction of the x's with the Q's are, in general, non-linear in X(t). For the 
special coarse grainings in which the variables are divided into a set distinguished by the 
coarse graining and a set ignored, we therefore have a general derivation of the form of the 
phenomenological equations of motion. We now discuss in more detail the implications of 
quantum-mechanical causality and quantum noise. 

(c) Quantum-Mechanical Causality Implies Classical Causality 

Feynman and Vernon [12] used path integral arguments to show that if £(t) is set to 
zero for t > £*, then W[X(r), £(r)] is independent of X(t) for t > £*. This result could be 
used to show that the force (F(t,X(r)]) is retarded, that is independent of the path X(t) 
for values of r greater than t. However the result also follows easily from the definitions 
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(fifty), (fiftytwo), and (fiftythree). In the Schrodinger picture 

(F(t,X(r)]) = Sp{u Ttt [X(T)] F(X(t)) U tfi [X(r)} p B 4 [X(r)}\ . (sixtyt 



wo) 



Using the cyclic property of the spur, the composition law, and unitarity of the evolution 
operators defined by (fortyone), we may write eq. (sixtytwo) in the form 

(F(t,X(r)}) = Sp{F(X(t)) U tfi [X(r)] p B u} fi [X{T)]} . (sixtythree) 

Since Utfi[X(r)} depends on X(r) only for < r < t, [cf. (fortyone)], this shows that 
(F(t,X(r)]) is retarded. 

The expression (one) for the decoherence functional incorporates a quantum-mechanical 
notion of causality. At one end of the histories, information about the specific closed sys- 
tem in the form of a density matrix must be supplied. At the other end is a unit matrix 
in the form of a 5-function representing a condition of indifference with respect to states 
at the end of the histories. This asymmetry between the two ends of the histories is the 
arrow of time in quantum mechanics. (See, e.g., [43, 2].) It is by convention that we call 
the extremity of the histories next to the density matrix p the past and consider time as 
increasing from it. 

To predict the future in quantum mechanics, we need the initial p and information 
about histories up to the present. We need no information about the future. That is 
quantum-mechanical causality. The retardation of (F(t,X(r)]) expressed by (sixtythree) 
shows that we need know the trajectory of the distinguished system only from the time 
of the initial p to the present to predict the system's future. That is classical causality. 
Eq. (sixtythree) thus shows that quantum mechanical causality implies classical causality. 

The origin and implications of the arrow of time in quantum mechanics may be use- 
fully discussed using a hypothetical generalization of the quantum mechanics of closed 
systems that employs both initial and final conditions [5, 3, 24]. In this generalization the 
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decoherence functional would be given by 

D(a'a)=N Sq / Sq Pf(q'f, qf) exp{i(S [q'{r)] - S [q{r)})/h\pi (q' , qo) , 

(sixty foura) 
where 

N~ X = dq dq' pf(q',q) Pi(q,q') ■ (sixty fourb) 

This generalization of quantum mechanics would permit the future and the past to be 
treated similarly. Arrows of time would arise in particular universes where p t is different 
from the time-reversed version of pj. In particular, usual quantum mechanics, represented 
by (one) and its associated causality and arrow of time arise for those universes, like ours, 
where pf oc / is a good representation of the final condition. In more general situations, 
with pf not proportional to /, there would be neither a notion of quantum mechanical 
causality nor a notion of "state of the system at a moment of time". Nor would the 
argument described above succeed in deriving classical causality. There would be advanced 
as well as retarded effects. 

(d) Quantum Noise 

In Section V we showed how the probabilities of decoherent sets of coarse-grained 
histories of the type under discussion could be thought of as the probabilities of a classical 
system in which the distinguished coordinates x a evolve from probabilistically distributed 
initial conditions according to an equation of motion in the presence of noise. In the case 
of a factored initial density matrix, initial conditions and noise are separately distributed 
and the system may be thought as obeying the Langevin equation (tenne), 

R(t) = £(t, X(r)} + C(t, X(r)] = , (sixty five) 

where the total force R(t) is distributed according to the generally non-positive distri- 
bution g[R(r),X(r)] constructed from the decoherence functional according to (tennd). 
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The analysis of this section provides an explicit form for £(t,X(r)] in (sixtyone) and a 
systematic approximation scheme for the spectrum of the noise. 

When decoherence is good enough that the restriction on the range of the ^-integration 
in (fiftyseven) arising from the coarse graining can be neglected, a systematic perturbation 
scheme for the approximate probabilities p(a) can be obtained by expanding the exponent 
in powers of £(£). In the leading approximation (sixty), an explicit expression for the 
distribution of the Langevin force is obtained, which is 

T T 



[det(.Kj/27r)]~2 e xp 



\ I dt J dt'c\t,X{T)]Kf Y {t,t'-X{T)}C{t,X{T) 



% 





(sixtysix) 



In this leading approximation, the noise is distributed with a positive Gaussian distribution 
function whose spectrum is fixed by the correlation function 

(jC(t,X(r)}jC(t',X(r)]) c = (h^K^t'; X(r)\ 

^^({AF(t,X(r)],AF(t',X(r)}}) . (sixtyseven) 

For the linear problem, that approximation is exact [cf. (twentyd)], with the further sim- 
plifying feature that the spectrum of Langevin force is independent of the path X(t). 

Higher order terms in the expansion of the exponent of (fiftyseven) may be regarded 
as providing corrections to the Gaussian noise. The general expression for the correlation 
functions is (tenk). Other corrections to the Gaussian noise arise from the coarse-graining 
restrictions on the range of ^-integration. 

VIII. Examples and Comparison with Classical Cases 

In this Section we specialize the general non-linear theory of the preceding Section to 
some particular cases considered by ourselves and other authors. We begin with the linear 
models described in Sections VI. 
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(a) Recovery of the Linear Theory 

In Section VI we derived the equations of motion for linear systems defined by an influ- 
ence phase (twelve) that was quadratic in the coordinates, x, distinguished by the coarse 
graining but containing no linear terms. An influence phase of this form will arise when 
the free action is quadratic in the coordinates distinguished by the coarse graining as in 
(eleven) when the interaction with the ignored coordinates Q is linear as in (thirteen), and 
when the density matrix has the special form (thirteena). To recover the linear equation of 
motion (nineteen) from the non-linear (fiftyeight), we need to evaluate (F(t, X(t)]) under 
these conditions. Eq. (sixtythree) gives a general expression for (F(t,X(r)]). We know 
from the general arguments described in Section VI that, when the interaction between dis- 
tinguished and ignored variables is linear as in (thirteen), the influence phase is quadratic 
in the x's and (F(t,X(r)]) is linear in X(r). It therefore suffices to evaluate (sixtythree) 
to linear order in a perturbation expansion in X(r); the higher orders must cancel. The 
result is 

(F(t,X(r)}} = (/(i)) + (i/ft) f dt'([f(t)J(t')]} X(t') . (sixtyetght) 

Here, the forces f(t) are the Heisenberg picture representatives of the homogeneous, linear 
functions of the ignored coordinates f(Q(t)) [cf. (fortynine)], and the expected values ( )o 
are computed using the time dependence of operators provided by Hamiltonian Hq of the 
ignored coordinates Q alone, neglecting their interaction with the x's [cf. (thirteenb)]. If 
we assume, as we did in Section VI, that (/(t))o = (this may always be achieved by 
a time-dependent shift in x(t)) then we recover both the form of the linear equation of 
motion (nineteen) and the expression for the additional force anticipated in (thirty). 

The derivation of the expression (thirteenb) for the spectrum of quantum noise in the 
linear problem is even more straightforward. From the form (fiftynine) we see that Kj(t, t) 
is equal to kj(t,t) of the linear case. Since we assumed {f(t)) vanished for the linear 
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problem [cf. (twentynine)], eq. (fiftynine) involves only the expected value of {f(t),f(t')}. 
However, general arguments for the linear problem show that when the interaction is of the 
linear form (thirteen) this expected value, proportional to kj(t,t), is independent of the 
x's. It may thus be evaluated in zeroth order in perturbation theory in x(t); that is, the 
time dependence of the operators in (fiftynine) is provided by the Hamiltonian Hq with no 
interactions between Q's and x's. This is the time dependence we denoted by a subscript 
zero in Section VI, and thus we recover (thirtyone). 

(b) Semi-Linear Systems 

Simple expressions for the spectrum of quantum noise and the kernel of the associated 
fluctuating force analogous to eqs (thirty) and (thirtyone) have been obtained in more 
general cases than strict linearity A particularly simple case occurs when the dependence 
of the total action on the ignored variables (the Q's) is restricted to be quadratic, but 
arbitrary dependence on the distinguished variables (the x's) is allowed in the potential 
energy and interaction terms. This case was discussed in purely classical situations by 
Zwanzig [25], and in particular field theory examples by Ryang and Saito [26], and more 
recently in a general survey of such problems by Brun[27], who applies our methods. 

Let us consider a specific class of problems. For simplicity, we assume a single distin- 
guished variable, x, to avoid matrix notation. Suppose that (i) the action of the ignored 
variables So[Q(t)] is a quadratic functional of the Q's; (ii) the interaction of the Q's and 
x's is of the form 

Suit [ x (t)i Q(t)\ — — / dt a(x(t))f(Q(t)) , (sixtynine) 



'0 

where / is homogeneous and linear in the Q's, but a is not necessarily a linear function of 
x; and (iii) the density matrix factors as in (eight) and has the special form (thirteena). 
Under these conditions the influence phase will be a functional that is at most quadratic 



We thank T. Brun for these references. 
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in a(x(t)), and its form is given by (twelve) with x(t) and x'(t) replaced by a(x(t)) and 
a(x'(t)) respectively, assuming, as before, that linear terms vanish. The equation of motion 
will then be of the form (fortysix), with the force arising from the interaction given by 

(F(t,X(r)]) = a'(X(t)) [ dt 1 k R {t,t') a (X(t')) , (seventy) 

Jo 

where a' = da/dX. Here, k R , as before in (thirty), is proportional to the expected value 
of the commutator of f(Q(t)) evolved in time without the interaction with the x's. Linear 
terms in the influence phase would give rise to a non-zero value of (/ (t))o an d an additional 
term in (seventy) of the form a (X(t))(f(t))Q. With the addition of this term, the result 
(seventy) is equivalent to that of Brun [27], who does not eliminate linear terms in the 
influence phase. The Langevin force C(t) is given by the formula 

C(t;X(r)} = a'(X(t))£(t) , (seventyone) 

where £(t) is a Gaussian random force with a spectrum given by (seventyone) and (twenty- 
five) . The Langevin equation describing both dissipation and fluctuations in this approxi- 
mation is then 



S(t) 



dSfreeiXi^ 1 '' 



^1 + a!{X(t)) [ dt'k R (t, t ! )a (X(t')) + £(t) 
JO 



. (seventya) 



SX(t) 
A fluctuation-dissipation theorem of the form given by (thirtysix) continues to hold. 

(c) Equations of Motion for Semi-linear Classical Coarse-Grained Systems 

The evolution of a classical Hamiltonian system is, of course, deterministic and essen- 
tially reversible when followed in all detail. In a coarse-grained description, however, such 
systems will, in general, approximately obey effective equations of motion including the 
effects of dissipative forces and with deviations from these equations produced by classical 
noise. The classical problem analogous to that considered in this paper would be to derive 
the Langevin equation for the motion of some followed variables x(t) that are interacting 
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with some ignored variables Q(t) whose initial conditions are probabilistically distributed 
according to some given rule. To our knowledge, this kind of problem has not been con- 
sidered classically for the non-linear situations discussed in Section VI. However, it has 
been worked out by R. Zwanzig [25] for the semi-linear systems treated in the preceding 
subsection. T. Brun [27] has shown that the Langevin equation deduced from quantum 
mechanics coincides with Zwanzig's result in the limit h — ► as it must. We briefly review 
Zwanzig's derivation and Brun's demonstration here. 

Again, we assume for simplicity a single distinguished variable x(t) interacting with 
an assembly of oscillators according to (sixtynine). The free action So[Q(t)] and the func- 
tion f(Q(t)) are given explicitly by (thirteenb) and (thirteenc) respectively. The classical 
equation of motion for the Q's following from (thirteenb) and (sixtynine) may be solved 
explicitly with the result 

1 rt 



Q A (t) 



oA 



H) 



dt' sin [w A (t - t')] C A a(X(t')) 



(seventyoneb) 



mu A j Q 

Here, I (t) is the solution of the free oscillator equations following from the action (thir- 
teenc) with the same initial position and momentum as Q it). This result for Q (t) may 
be substituted into the classical equation of motion for X(t). The result is an equation of 
motion for X of the form 

-t 



^S ircc [X(r)} , 

6X(t) +a{X{t >> 



dt'k R {t,t')a{X{t')) + £(t) 







(seventyonec) 



where 



[seventyoned) 



£(t) = J2 c Ant) 

A 
and we have made use of the fact that the sum of retarded Green's functions from (seven- 
tyoneb) enters (seventyonec) as the combination we called kji(t,t') in (thirteene a). 

If the initial conditions for the Q's are probabilistically distributed, then the motion 
X(t) will be probabilistically distributed as well. The initial conditions of the Q's are the 
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initial conditions of the free oscillator motions £ it). In (seventyonec), the only way these 
initial conditions enter is through the function £(t) defined in (seventyoned). Equation 
(seventyonec) may, therefore, be interpreted as a Langevin equation with a stochastic force 
a'(X(i)£(i). The time dependence of this force is known because the I it) in (seventyoned) 
satisfy the harmonic oscillator equations of motion. When the values of the Q's and their 
conjugate momenta are distributed thermally, the correlation functions of this noise are 
easily calculated. The classical phase space distribution analogous to the thermal bath 
used in [11] and [12] is 

p d (n , Qo) oc exp [-Hq(Uq, Qo)/kT B ] , (seventyonee) 

where Hq is the classical Hamiltonian corresponding to Sq. (A slightly different distribution 
was assumed by Zwanzig in [25] with slightly different results.) The result is Gaussian noise 
with (£(t)) c = and 

(£(t)£(t')) c = tfV -^cos \oj A (t-t')] = -kf{t,t ! ) . (seventyonef) 

A muJ A ~ 2 

Eq (seventyonec) and (seventyonef) are essentially the results of Zwanzig [25]. We note 
that the spectral weights of the spectrum of the fluctuations (seventyonef) and of the kernel 
of the dissipative force term (thirteene) as defined by (thirtyfive) are related by 

kf(u>) = (kT^/u)kji(u) . (seventy oneg) 

This is the classical fluctuation-dissipation theorem. 

As discussed by Brun, the Langevin equation of the quantum-mechanical problem 
(seventya) coincides in form with (seventyonec) derived classically by Zwanzig [25]. The 
only difference is that the noise spectrum hkj(t, t) is given by (thirteene b) in the quantum- 
mechanical case and by (seventyonef) in the classical one. Indeed, the two expressions 
coincide in the classical limit as they must. The classical noise is entirely thermal. When 
quantum mechanics is taken into account there is quantum noise as well. 
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IX. More General Coarse Grainings 

The coarse grainings discussed in the previous sections are limited to those that dis- 
tinguish a fixed subset of the coordinates q l . Coarse grainings that realistically describe a 
quasiclassical domain are not of this simple type. As we have discussed elsewhere [2] , it is 
likely that a quasiclassical domain will be described by, among other things, coarse graining 
with respect to ranges of values of the averages of densities of conserved or approximately 
conserved quantities over suitably small volumes. Examples are the densities of energy, 
momentum, charge, current, nuclear species, etc. Together with field averages, these are 
the "hydro dynamic" variables that enter into the differential equations of classical physics. 
Sufficiently large volumes would give these variables enough "inertia" to enable them to 
resist the deviations from predictability caused by the interactions that effect decoherence, 
as we have described in the earlier sections. 

The coarse grainings discussed in this paper must be generalized in two ways to discuss 
such variables. They must be generalized to allow the original fine-grained description to 
involve momenta as well as coordinates. They must also be generalized to permit coarse 
grainings by ranges of values of averaged densities, These correspond to no particular fixed 
subset of coordinates. In this section we introduce the machinery necessary to consider 
such coarse grainings, although we do not carry out an analysis of the circumstances in 
which they decohere or behave quasiclassically. 

(a) Phase-Space Coarse Grainings 

A fairly general class of coarse grained histories may be obtained by considering parti- 
tions of the coordinates q 1 at some times and their conjugate momenta 7Tj at others. More 
specifically, we consider partitions of configuration space into an exhaustive set of exclusive 
regions {A^ } at some of the times t^ and partitions by partitions of momentum space 
into an exhaustive set of regions of momentum space {A*j, } at other of the times t^. In 

74 



this Section, concerned only with this type of coarse graining, we will reserve the notation 
{Pa k {tk)} for projections of configuration space alternatives and introduce the notation 
{Pa k {tk)} for a set of orthogonal projections onto a set of exclusive momentum space re- 
gions. The sets of histories we are considering consists of sequences of sets of either P's 
or P's at times t\, ■ ■ ■ , t n . The decoherence functional for such a set of histories is given 
generally by (fourfive). We now show that it has a sum-over-histories representation by 
path integrals in phase-space. 

Utilizing complete sets of coordinate eigenstates, we may write the decoherence func- 
tional (fourfive) in the form 

D(a',a) = dq' f dq f / dq / dq 5(q' f - q f) (q f T\C a i\q'$) p (%,qo) (qoO\ci\q f T) . 

(seventy five) 
The matrix elements of the C a may be written as the compositions of sequences of propa- 
gators between definite coordinate or momentum eigenstates. For example, if P^^l) is a 
projection onto a momentum region A^^ti) and P& 2 (t2) onto a coordinate region A^fo) 
we can write 

(qfT\C a \q Q 0) = (Q / T|P^ 2 (t 2 )Pa 1 (ti)|Q 0) (seventy six) 

in the Schrodinger picture as 

/ dq 2 I dni (qfT\q 2 t 2 ) (<?2*2Kl*l) (^l^lko ) • (seventy seven) 

The propagators in (seventy seven) may be represented as phase-space path integrals if the 
Hamiltonian H(iT,q) associated with the action £[(/(£)] is of a suitably simple form [28, 
29]. In particular, they can be so represented if H(ir, q) is a sum of a function of the 7r's 
and a function of the g's. For example, 

(q t \irt) = I 5fj, ex.p{iS[n(t),q(t)]/h} , (seventyeight) 
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where S[tt, q] is the canonical form of the action: 

ft" 
S[TT,q]— / dt[wi(t)q (t) — H(7r(t), q(t))] . (seventy nine) 

Jt> 

The integral is over phase-space paths between t and t weighted by the invariant Liouville 
measure and restricted by the conditions that they intersect n at t and q at t . The 
details of these integrations are spelled out in the Appendix. 

By inserting (seventynine) into (seventyeight) and (seventyeight) into (seventysix) , one 
arrives at a sum-over-histories form for the decoherence functional for phase-space coarse 
grainings 

D(a, a)= ISfi' / 6fi S(y - qj) exp|i(5[7T / (r), q'(r)] - S[w(r), q(T)])/hjp(q , %) . 

(eighty) 
The integral is over phase-space paths restricted by the coarse graining. For example, the 
integral over n(t) and q(t) is over the phase-space paths that thread the intervals in either 
coordinate or momentum space corresponding to the history a. There is no integration 
over the initial momenta 7Tg and ttq but unrestricted integrations over the final momenta 
jv'r and TTf. Again, the details of this and the measure are in the Appendix. 

In usual cases where the Hamiltonian in (seventynine) is quadratic in all momenta, the 
momentum dependence in the integrand in (eighty) is that of a Gaussian. For those coarse 
grainings that restrict only the coordinates q l and ignore the momenta 71^, the integrals 
over momenta may be carried out explicitly. The result is the Lagrangian path integral for 
the decoherence functional (one). Indeed, it is by this route that the measure in that path 
integral is usually derived from the canonical, Liouville, u dpdq/(2Trh)" measure on paths 
in phase-space. 

There is no obstacle to letting the time of a coarse graining by momenta coincide with 
that of a coarse graining by coordinates. Even in the quantum mechanics of measured 
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subsystems it is possible to consider a measurement of position followed after an arbitrar- 
ily short time interval by a measurement of momentum. Care must be taken, however, to 
specify the order of the coarse grainings when two such times coincide. Since the corre- 
sponding operators do not commute, a projection on a range of momentum at one time 
followed immediately by a projection on a range of position defines a different history from 
one with the operators in the opposite time order. 

With this machinery in hand, we may now consider phase-space coarse grainings, 
analogous to those of Section III, in which the phase-space coordinates (7^, q l ) are divided 
into cannonically conjugate pairs (p a , x a ) that are distinguished by the coarse graining 
while the remaining pairs (II^, Q ) are ignored. A simple, interesting, and important 
class of models is obtained by assuming that the action decomposes according to (3.2) 
with Sfj-ee and Sq quadratic in time derivatives and the interaction independent of all time 
derivitives. (A slight generalization would then be needed to deal with a system of particles 
interacting electromagnetically.) The decomposition of the Hamiltonian corresponding to 
these assumptions is 

H(n, q) = H ircc (p, x) + H (U, Q) + H int (x, Q) , ((eightyone)) 

where -fff ree and Hq have a quadratic momentum dependence. The simplifying consequence 
of these assumptions is that the Gaussian integrals over the momenta II^ in (eighty) may 
all be carried out explicitly. The remaining integrals over the Q have the same form as 
they do in the Langrangian path integral and may be summarized by a single influence 
phase W[x'(t),x(t)] defined by (three). The decoherence functional for phase-space coarse 
grainings may then be written 



D{a ,a) = / d/j' / dfi S(xf — x'A 
J a' J a 



xexp|i(^ free Ip(t),x(t)] -S bcc [p(T),x(T)] + W [x'(t),x(t)] )/hjp(x' ,xo) , 

(eightytwo) 
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with p defined by (four). 

Eq. (eightytwo) shows that coarse grainings in which a fixed set of coordinates or their 
conjugate momenta are followed at a sequence of times while all others are ignored may 
be studied by path integral techniques. The form of (seventynine) , however, allows an 
immediate and important distinction to be drawn between coordinate coarse grainings 
and momentum coarse grainings. Coordinate alternatives will decohere if W has a positive 
imaginary part that becomes large as x (t) and x(t) are increasingly distinct. However, 
there is no corresponding mechanism leading to the decoherence of momentum alternatives 
for this class of models. 

In the case of the linear models discussed in Section VI, this conclusion may be made 
more precise by following a few steps that led to the derivation of the equations of motion. 
Introduce variables for the momenta analogous to those for the coordinates in (tend) 

zu(t) — p (t) — pit) , (eightythreea) 

I 

P(t) = - [p(t)+p(t)] . (eightythreeb) 

The exponent in (eighty) may now be reexpressed in terms of the variables of (tend) and 
(eighty) and, after a few integrations by parts, put in the form 



%cc [p'(t),x'(t)] + S bcc [p(r), x(t)] + W [x'(t),x(t 
-&0 + £ dt{w\t) [X{t) - J£L] + e + (t) [-P(t) + ^y + £ dt'k R {t, t')X{t')] } 



»T r T 



+ - I dt f dt'^^k^t')^') . (eighty four) 

4 JO JO 



The terms in square brackets on the right hand side of (eightyfour) are Hamilton's equations 
of motion augmented by terms describing the additional forces arising from the interaction 
of the (p, x) subsystem with the rest. However, we cannot conclude that these equations 
of motion are valid for phase-space coarse grainings. The last term in (eightyfour) makes 
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the integrand of (eightytwo) small when £ 7^ and thereby enforces the decoherence 
of coordinate alternatives. However, the absolute value of that integrand is uniformly 
distributed in w. Unless the integration over the other variables makes (eightyfour) small 
for w 7^ there will be no decoherence of momentum alternatives. In cases when 
momentum alternatives cannot even be assigned probabilities there is a fortiori no issue 
of whether the probability is high for their correlation in time by an effective equation of 
motion. The origin of this distinction between coordinates and momenta in this model 
is, as has been remarked by many authors, that the interaction Hamiltonian is local in 
coordinates but not in momentum. 

This analysis of both coordinate and momentum coarse grainings stresses an impor- 
tant if familiar point. Two coarse grainings that would be essentially indistinguishable 
classically may have very different properties quantum-mechanically. The present exam- 
ple illustrates this with coarse grainings by momentum, p(t), and by the difference in 
coordinates at two nearby times M[x(t + e) — x(t)]/e. As we have seen, coarse grain- 
ings by the momentum even at one time do not decohere, while coarse grainings by the 
coordinates at times separated by intervals longer than the decoherence time-scale set 
by kj [cf. (fourthirty seven)] do decohere. Thus, these two types of coarse graining have 
essentially different properties with respect to decoherence even though they would be in- 
distinguishable on the basis of classical physics when the decoherence time-scale is much 
shorter than characteristic dynamical time-scales of the classical motion. Evidently, con- 
siderable care is required in identifying the variables through which quasiclassical behavior 
is to be defined in quantum mechanics. At the present moment it is, therefore, an open 
question for the purpose of defining classicality whether a quantity like momentum density 



Such cancellation leading to the decoherence of momentum alternatives occurs, for exam- 
ple, in the case of a free, non-relativistic particle. The conservation of momentum means 
that projections onto ranges of momentum at different times commute and a history com- 
posed exclusively of such projections will automatically decohere. We thank J. Halliwell 
for pointing this out to us. 
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is best defined in terms of time derivatives of fields or in terms of differences of fields over 
a macroscopically negligible interval. 

(b) Densities 

The most important generalization of the coarse grainings studied in this paper is to 
cases where the variables that are distinguished are not limited to a fixed set of funda- 
mental coordinates. To coarse-grain by the value of the baryon number in a small spatial 
volume, for example, is not the same as following some particular subset of the fields of 
a fundamental field theory. We shall discuss how, in principle, more general and realistic 
coarse grainings can be treated by techniques analogous to those used to discuss the special 
cases of this paper. Our considerations are essentially formal and we have not pushed the 
analysis far enough to derive equations of motion in these realistic cases. Our discussion, 
however, indicates a route by which that might be accomplished. 

A very general and useful class of coarse grainings is obtained by partitioning the fine- 
grained histories according to the values of functionals of them. To illustrate an interesting 
case in a manageable notation, we consider a field theory with a single charged scalar field 
<f)(x,t). A set of functionals leading to coarse grainings relevant for the present discussion 
consists of the values of the charge density at time t averaged over a small spatial volume 
v at different spatial points x: 



J°(x,t) = - [ d 3 £ j°{x + i,t) (eighty five) 

v J v 



where j (x,t) is the charge density expressed in terms of (j)(x,t). The fine-grained his- 
tories, <f)(x,t), may be partitioned into exhaustive and exclusive classes by the values of 
J (x, t). (Partition by ranges of values of these averaged densities would be a further 
coarse graining.) Consider a particular value v(x, t). (The conventional use of p for charge 
density is precluded by its use here for the density matrix!). The coarse-grained history 
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corresponding to v(x,t) consists of all those cf>(x,t) for which the integral (eightyfive) has 
this value. 

The decoherence functional for a pair of histories coarse-grained by particular values 
of the charge density is 

D[v'(x,t),v(x,t)] = [ 5</>' fs<f>S[<f>'(x,T)-<l/(x,T)] 



x5 



J ui (x,t)-v'(x,t) exp\-(S[(f>'(x,t)]-S[(l)(x,t)])\6 J U (£, t) - v(x, t) p[<j/(x,0),<l>(x,0)\ 



rO/ 



(eightysix) 
The first 5-functional in (eightysix) enforces the coincidence of the histories at the final 
time T as in (seventyfive) . The other two 5-functionals restrict the fields in the integral to 
the coarse-grained histories labeled by v'(x, t) and v(x, t). Of course, making precise sense 
of a formal expression like (eightysix) raises many mathematical issues that we shall not 
pursue here. 

The decoherence functional D[v (x, t), v(x, t)] is the generalization of D[x (t), x(t)] con- 
sidered in the simple models of the earlier sections of this paper. In those cases a fixed set 
of variables describing the fine-grained histories of a single oscillator were distinguished 
by the coarse graining while the variables describing the other oscillators were ignored. In 
the case of a coarse graining by the average value of a density over a volume, one could 
loosely say that the "variables" describing the different field configurations internal to the 
volume and consistent with a given average value are ignored. Hence there is no subset 
of the fundamental fields or linear transformation of them that describes these "internal 
variables". That is why D[v'(x,t),v(x,t)] is a generalization of D[x'(t),x(t)]. 

The form of D[v (x,t),v(x, t)] may be brought in closer analogy with (five) by repre- 
senting some of the 5-functionals in (eightysix) as exponential integrals 



8[J (x,t) — v(x,t)] = J Syexp \i j d x y(x,t) [J (x,t) —v(x,t)] j- . (eightyseven) 
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Then, eq. (eightysix) becomes: 

D[v'{x,t),v(x,t)] = I Sy' f 5y f 5</>' f5(j)5[(j)\x,T)-(j)(x,T)]exp(Us[(j)'(x,t)]-S[(j)(x,t)} 

+ I d xy(x,t) J (x, t) — v (x, t) — I d x y(x,t)[J (x, t)— v(x, t)] J >p [4> (x, 0),(f>(x, 0)] . 

(eightyeight) 
The integral in (eightyeight) is of a familiar form in a field theory with sources and should 
be accessible to standard approximation techniques. Of course, to study the putative 
"hydro dynamic" variables of a quasiclassical domain, coarse graining by the densities of 
other approximately conserved quantities such as energy, momentum, baryon number, 
as well as the averages of long-range fields, must be considered. Once the decoherence 
functional is calculated from expressions like (eightyeight), the circumstances in which 
such coarse-grained sets of histories decohere can be investigated and their equations of 
motion derived by the methods of Section VII. An important question will be the closure 
of any such set of equations of motion. However, we have not yet progressed beyond this 
formal sketch of the route to a derivation of the classical equations of motion for realistic 
quantum systems. 

X. Conclusions and Program 

Instead of merely summarizing the results obtained in this article and the future di- 
rections of research sketched in the last Section, let us relate them to the program (of 
understanding quantum mechanics and its relation to quasiclassical experience) in which 
we have been engaged, and which involves many elements elucidated by other authors over 
the last thirty-five years. We shall briefly review that program and where the present work 
fits into it. 

We start with the quantum mechanics of a closed system representing the universe 
(with a deliberately simplified treatment of the complications caused by quantum grav- 
ity). There is a dynamical theory of all the elementary particles and their interactions, 
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which we take to be described by an action or a Lagrangian or a Hamiltonian H, and 
also a Heisenberg density matrix (or initial density matrix) p. The question is how the 
quasiclassical domain of familiar experience comes to be an emergent feature of the system 
characterized by H and p. 

A "measurement situation" can often be characterized as one in which some variable 
comes into strong correlation with the quasiclassical domain. The quasiclassical domain 
also permits, for certain probability tracks, certain spatial regions, and certain epochs 
of time, enough classical predictability for the evolution of complex adaptive systems 
that learn, observe, and record, and also utilize, in some approximation, the probabil- 
ities assigned by quantum mechanics, on the basis of p and H, to different alternative 
coarse-grained histories. (Presumably, it is the quasiclassical correlations, representing 
near certainties, that are easiest to utilize.) Observations by such complex adaptive sys- 
tems, which we call IGUSes (information gathering and utilizing systems) when they are 
functioning as observers, are then considered to be actual measurements. 

A quasiclassical domain is defined by decoherence, a measure of classical predictability, 
and some sort of maximality (such as what we have called "fullness" [3]). The first re- 
quirement is decoherence, that is, enough coarse graining of the alternative histories of the 
universe so that there is, exactly or to a very good approximation, no interference between 
the alternative coarse-grained histories, as measured by the decoherence functional. 

Decoherence requires coarse graining that goes far beyond the modest requirement im- 
posed at each instant of time by the uncertainty principle. The mechanism of decoherence 
involves the loss of phase information as a result of the coarse graining and is associ- 
ated with noise that inextricably combines quantum fluctuations and classical statistical 
fluctuations, both necessary for the decoherence. 

All of that is most obvious, of course, in the instructive, although oversimplified models 
in which the coarse grainings include an average over some fixed set of variables (the ignored 
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ones ) while following others (the distinguished ones) more or less coarsely. The ignored 
variables, through their interactions with the distinguished ones, carry away the phases 
and are responsible for the fluctuations. 

The fluctuations, of course, cause departures from any effective (or phenomenological) 
classical equations of motion for the distinguished variables. But a high degree of deco- 
herence requires very large fluctuations, which threaten to produce very great departures 
from classical predictability. Thus, for the coarse graining to yield a quasiclassical domain, 
it is essential that the distinguished variables carry very high inertia so as to resist most 
of the large fluctuations and follow the effective equations of motion with only small devi- 
ations over long stretches of time and with only occasional large ones. That high inertia 
is achieved by even much coarser graining than was required for the decoherence itself. 

The requirement of fullness, a kind of maximality, was discussed in [33]. For the 
coarse graining defining a quasiclassical domain to be an emergent feature of the universe 
characterized by H and p rather than an artifact chosen by some IGUS, it should be 
as refined a description of the universe as possible consistent with the requirements of 
decoherence and quasiclassicality. In [33] we proposed, for the case of perfect decoherence, 
the notion of a "full" set of alternative coarse-grained histories to capture this idea of 
maximal refinement. Any "full" set belongs to an equivalence class defined by a basis 
in Hilbert space, provided the density matrix p corresponds to a pure state. (If it does 
not, and especially if it has a great many eigenstates with non-zero eigenvalues, then 
the condition of "strong decoherence" discussed in [33] is too strong and the associated 
discussion of maximality requires modification.) 

We have posed the question as to whether there could be various kinds of essentially 
inequivalent quasiclassical domains or whether any quasiclassical domain is more or less 
equivalent to any other. The former case poses some challenging intellectual puzzles, 
especially if we imagine IGUSes evolving in relation to each of the essentially inequivalent 
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quasiclassical domains. 

There is, of course, an indication from our everyday experience of some of the features 
of a particular quasiclassical domain. It involves distinguished variables that are more or 
less hydrodynamic in character; they are integrals or averages over small regions of space 
of conserved or nearly conserved densities and of fields coupled to those densities. The 
regions must be large enough to produce sufficient inertia to resist most of the fluctuations 
associated with the coarse graining and small enough to implement the requirements of 
"fullness" . It is clear, however, that the distinguished variables cannot be defined once and 
for all, but depend on history. For example, the suitable hydrodynamic variables under the 
conditions that prevailed before the condensation of the solar system involved much bigger 
volumes than those that were suitable inside the planets after they were formed. That is 
why the models in which the ignored variables are fixed for all time are only instructive 
examples and not general enough for the realistic case. 

Our program thus aims at describing a quasiclassical domain with history-dependent 
distinguished variables resembling hydrodynamic ones and obeying effective classical equa- 
tions of motion apart from small fluctuations and occasional large ones, some of which 
result in the need for redefinition (for later times) of the distinguished variables. 

In this article we have, for the most part, confined ourselves to a model in which 
the distinguished variables are separated once and for all from the ignored ones and are 
also, unlike the hydrodynamic ones, coordinate variables of the fundamental theory, like 
modes of a scalar field, with kinetic energy bilinear in the time derivatives and the rest of 
the energy not involving the time derivatives. We have also restricted attention to initial 
density matrices that factor into the product of a density matrix for the distinguished 
variables and one for the ignored ones. With those simplifying assumptions, we identify 
the effective classical equations of motion and the Gaussian part of the noise that disturbs 
them. 
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We divided into two parts the quantum-mechanical process of prediction for further 
coarse grainings of the histories of the distinguished system. First, there is the calculus 
of amplitudes for histories fine-grained in the distinguished variables and of the bilinear 
combinations of amplitudes that define the decoherence functional. Second, there are the 
requirements of decoherence for deriving consistent probabilities from these amplitudes for 
histories that are further coarse-grained. 

We assumed that the requirements of decoherence were satisfied (although we discussed 
mechanisms by which this happens) and examined when the probabilities for histories 
favored classical predictability. 

We started with a general formulation that applies to fully non-linear systems. By 
performing a functional Fourier transform on the decoherence functional with respect to the 
difference in distinguished coordinates on the left and right, we introduced a distribution 
function for the "total force" (including the inertial term) that acts on the distinguished 
system along its history. Like the analogous Wigner distribution at a moment of time, it 
is not generally positive. 

Utilizing the distribution of the "total force" , we were able, even in non-linear cases, to 
represent the probabilities of a decoherent set of coarse-grained histories as the probabilities 
of a classical dynamical system governed by a Langevin equation incorporating history 
dependent noise. The initial conditions of this classical system are distributed according 
to the Wigner distribution and the noise according to the distribution of "total force" 
mentioned above. Both distributions are generally non-positive, which distinguishes them 
from classical ones, although the resulting probabilities for decohering sets of histories 
are, of course, positive. When the noise is small compared with the inertial term in the 
equation of motion, the coarse-grained histories are classically deterministic. 

As suggested by the simple physical picture of phases being carried away by the in- 
teraction of the distinguished subsystem with the ignored variables, there is a connection 
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between decoherence and noise. This was exhibited by writing the decoherence functional 
for the distinguished variables as an integral of an average density matrix times an expo- 
nential, and expanding the argument of the exponential in powers of the generalized vector 
£(£) that measures the difference between the distinguished variables as functions of time 
on the left- and right-hand sides of the decoherence functional. After a partial integration, 
the linear term has the form i/% j £}{t)E(t)dt and the equation £(t) = is precisely the 
effective classical equation of motion for X(t), which is the average of the distinguished 
coordinate variables on the left and on the right. The second term in the expansion of 
the argument of the exponential is —l/A%jj^{t)Kj(t,t')^{t')dtdt', where Kj is the pos- 
itive kernel describing the Gaussian noise, and can be thought of as the self-correlation 
(depending in general on X(t)) of a Langevin force £(t) added to the equation of motion. 

There are, in general higher terms in the expansion, corresponding to the fact that 
coarse-grained quantum mechanics is not exactly equivalent to effective classical equa- 
tions of motion accompanied by Gaussian noise. However, we are interested in the case 
of approximate decoherence of the coarse-grained histories, meaning that £(t) is mostly 
confined to very small values. When Kj is very large, this is achieved, provided that the 
higher terms are not somehow still more important. However, large Kj means large noise 
and that is the connection between decoherence and noise. We, therefore, also impose the 
requirement of a very large inertia matrix M, so that the huge noise is mostly resisted 
and the effective classical equations of motion are followed with high probability by the 
distinguished variables. 

For pedagogical reasons, we started by reviewing the completely linear case, studied 
extensively by Feynman and Vernon, Caldeira and Leggett, Unruh and Zurek, etc. There 
the noise correlation kj(t,t) is just a numerical function and the effective equation of 
motion for X(t) is linear, with a frictional force J kji(t, t)X(t)dt that is in general nonlocal 
in time. When the effective density matrix is diagonal in the energy of the ignored variables, 
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then the numerical functions kj and kji depend on t — t only and are related to each other 
by the fluctuation-dissipation theorem. In the Fokker-Planck limit, kj is proportional to 
S(t — t') and kji to S'(t — t') and the coefficients are related by a very simple fluctuation- 
dissipation relation. 

We then generalized to the nonlinear case discussed above. That is the principal content 
of this paper, and the work can be applied to many problems, as in [1] and [27]. 

We also treated more general coarse grainings, in which both generalized coordinates 
and generalized momenta are utilized and what had been a Lagrangian formulation turns 
into what is essentially a Hamiltonian formulation. 

Finally, we pointed out the desirability of removing the remaining two unrealistic 
conditions, that the distinguished variables have their time derivatives occurring in the 
action only in a bilinear kinetic energy term, and that the distinguished variables be 
defined independent of history. We are studying the mathematical problems posed by 
these generalizations. If the first difficulty can be overcome, then one should be able to 
treat the effective classical equations of motion and the Langevin force for distinguished 
variables that include the hydrodynamic quantities of familiar quasiclassical experience. 

If the second and greater difficulty can be overcome, and history dependence introduced 
into the coarse graining, then we may begin to tackle the deep problem of introducing 
individuality into quantum mechanics. Actual alternative histories deal, of course, in large 
part with individual objects such as our galaxy, the sun, the earth, biological organisms on 
the earth, and so forth. Yet discussions of quantum mechanics up to now have typically 
treated such individual objects only as external systems, labeled as "observers" and "pieces 
of apparatus" . If history dependence can be properly introduced into the explicit treatment 
of quantum mechanics, then we may be able to handle individuality with the care that it 
deserves. 



Appendix — Path Integrals 

In the body of this paper we have formally manipulated both configuration space and 
phase space path integrals. We define those integrals more carefully in this Appendix and 
use the definition to show that the formal manipulations we have used are legitimate. 

(a) Phase Space Integrals 

We begin by deriving an explicit expression for the most general path integral oc- 
curring in the decoherence functional for sets of histories coarse-grained both by regions 
of configuration space and regions of momentum space, as discussed in Section IX. To 
keep the notation manageable, we shall assume for the moment that we are considering a 
one-dimensional system with coordinate q, momentum n, and a Hamiltonian of the form 

H(n, q) = — + V(q, t) . (etghtymne) 

The generalization to larger-dimensional configuration spaces is obvious. A set of histories 
consisting of chains of just two projections, one on momentum space region {A Q , 1 } at 
time t\ and the other on configuration space region {A a2 } at time £2, will suffice to 
illustrate the construction in more general situations. The important matrix elements for 
the construction of the decoherence functional are of the form (seventysix). They may be 
expressed in the Schrodinger picture as 

(q f T\C a \q 0) = {q f T\P^{t 2 )P^{ti)\q^) 

= (qfle-^^-^^P^e-^^-^^P^e-^/^qo) . (ninety) 

Here we have written P& 2 and P^ to recall explicitly that we are dealing with projections 
onto a configuration space region A a2 and on a momentum space region A ai respectively. 

We now cast the right-hand side of (ninety) into a phase space path integral form. 
We shall be brief because the construction is a standard one [28, 29]. We divide the 
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interval [0, T] up into N equally spaced time slices to = 0, t, t, • ■ • , £jv-1> tjq = T with an 
interval e = T/N between them. We assume that the times t\ and £2 coincide with two of 
these slices for a sequence of iV's tending to infinity. Let K\ and K^ be the labels of the 
slices corresponding to t\ and t<i respectively, understanding that these are functions of N. 
Write the propagators in (ninety) as the product of an appropriate number of factors of 
ex.p{—iHe/h). Between these factors, on each time slice except the first, insert a resolution 
of the identity of the form 



/ dn k / dq k \q k >< q k \n k >< n k \ = I 



(ninetyone) 



The result is the following expression for the matrix element in (ninety): 

f N 

(q f T\C a \q 0) = / Y[dn k dq k (q f \q N )(q N \n N )(n N \e^ He / h \q N _i) 

J k=l 
x (qN-l\*N-l)(*N-l\e- iHe/h \qN-2) • • • (qK a \P&\*K a )(*K 2 \e- iHe/h \qK a -l) • • ' 

x (QK.lPaMK^^K^^'^lqK.-l) ■ ■ ■ (qikl)(^lK lHe/h \qo) ■ (ninetytwo) 



Now note the following relations: 

(q\ w ) = e i7Tq / n /(2>irh)5 , {ninety three) 

(q\Pl\n) = e a MJ* q,n /&**)* . {ninety four) 

(q\P^\ir) = e ai {n)e m( i/ h /{27rh) 1 2 , {ninety five) 

where e a {x) is unity for x in the interval a and zero outside it. (The symbol n stands for 
momentum unless it occurs in the combination 2irh.) Further, to first order in small e, 

(ir\e- iHe / h \q) w e- lH ^^ e l h e~ l7:q l h /{2ixh) l 2 {ninetysix) 

where H(jr, q) is the function given by (eightynine) . Inserting these relations in (ninetytwo), 
noting that there is a 5- function that identifies qjy and qt,we find the following expression, 
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which is exact as N — ► oo, 



mCam =^j^( N fi^\ 



e a 2 {(lK 2 ) e a 1 {^K 1 ) 



X 



exp uE 



N 



e 7r 



.7=1 



Qj ~ Qj-1 



H (TTj,qj-i,tj_i) 



{ninety seven) 



This is the definition of the phase space path integral that we have written in Section IX 



as 



(qfT\C a \qoO) = / SnSqexp 



h 



T 



dt 



n(t)q(t)-H(ir(t),q(t),t) 



(ninety eight) 



It is an integral over phase space paths in the class specified by the coarse graining, 
that is, over paths which pass through the momentum space region A ai at time t\ and 
configuration space region A a2 at time t<i- 

(b) Path Integrals for the Decoherence Functional 

When the coarse graining is only by regions of configuration space, and there is no 
coarse graining by momentum space, then the Gaussian integrals over the 7r&, k = 1, • • • , N 
may be carried out in expressions like (ninetyseven) leading to Lagrangian path integrals 
for the matrix elements (qtT\C a \qoO) corresponding to individual histories. For example, 
a history defined by a sequence of (/-intervals A ai , A a2 , • • • , A an at times t\, • • • , t n would 
have 



(qfT\C a \qoO) = lim 

J N^oo 



7V-1 / M \\ 

n {^k)^ 



k=\ 



ea n (QK n )---eai(QK 1 ] 



x exp 



h 



N-l 

E 

J=0 



l - M ( q i +l ~ q i 



V(q v t 



3i v 3> 



(ninetynine) 



Here, K^ is the label of the slice corresponding to t{ and, as before, qjy — qt. 
Eq. (ninetynine) is the sum-over-paths usually written 



(q f T\C a \q 0) 



[qoaqf] 
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5q exp(iS[q(t)]/h) 



(hundred) 



T 1 

-Mq z — V(q,t) . (hundredone) 



where S[q(t)] is the Lagrangian form of the action corresponding to the Hamiltonian 
(eightynine) . 

rp 

S[q(t)] = [ dt 
JO 

In an even more compact notation, where, as in (one), there is an integration over go an d 
qt, we have denoted the restrictions on the range of integration arising from the coarse 
graining by an unadorned subscript a on the integral sign. 

The path integral for the decoherence functional consists of two multiple integrals 
like (ninetynine) over two polygonal paths {q^} and {<?[,} with additional integrals over 
their initial and final endpoints weighted by the initial density matrix and final 5-function 
respectively. When, as discussed in Section III, the action is of a suitable form and the 
coarse graining constrains only a fixed subset x a of the variables of configuration space, the 
integrals over the remaining {Qf} may be carried out yielding a path integral involving 
an influence phase (five). To better understand how formal manipulations are carried out 
on that path integral we shall now write out the explicit time-slicing implementation of 
it following from (ninetynine). To keep the notation manageable we shall consider the 
case where the coarse graining refers to a single coordinate x and limit attention to the 
explicitly linear problem discussed in Section VI having a quadratic influence phase given 
by (twelve). The generalizations of this case should be obvious. The integral (five) for the 
decoherence functional is then explicitly: 

D(a , a) = lim / ( — | I dxudxi, 

1 ' N^ooJ L\2irieh) £A k k 

x5(x' N — xpj) E a / a (x f j, Xj) exp{iA(x'j, Xj)/Ti)\ p(x'q, xq) (hundredtwo) 

where the functions E a > a and A are defined as follows: The function E a i a is 

n 

Ea'a{ x 'ji x j) — J [ e a'X x> K i ) e o'i( x K i ) (hundr edthr ee) 

i=\ 
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and enforces the constraints of the coarse graining. The exponent, A, is the discrete form 
of (eleven) and (twelve), viz. 



2"-J 



2"^ 



N J 



H — >^ >^ e (x„- — Xj) [/c^(j, ^)(a^ + xg) + ikj(j, £)(xg — xg)\ . {hundredf our) 

j=ie=i 

Here, kji(j, £) and &/(,?', £) are the real and imaginary parts of the function k(t', t) evaluated 
on the discrete time slices. In passing from (ninetynine) to (hundredthree) we have made 
use of the fact that in the limit N — *■ oo it makes no difference in the Lagrangian path 
integral whether the integral J V(x(t))dt between time slices is approximated using the 
value of V at the start or end of the interval. The form (hundredfour) is slightly more 
convenient for what follows. 

We now change variables in the multiple integral (ninetynine) from {x'.^x^} to the 
discrete versions of (fifteen) 



^k — x k~ x ki ^k — n( x k + x k) • 



(hundredf ive) 



The Jacobian is unity on each slice so that (hundredtwo) becomes 

[\f M \ N ~ l N 

D{a '> a) = ^oo J [ \2m) n dxA 



X( 



mN)E a 'a (Xj + Co/2, Xj - Co/2) exp [iA(X jt ^)/n)] p(X Q + Co/2, X Q - £ /2) . 

(hundredsix) 

After a little algebra the exponent may be written 

N 



A(X V ^) = -^M 



+E^ 



j=i 



Xj+i — 2Xj — Xj_i 



+ KX 



N N 



+1 EE <?tj[*>RU>*)Xt + iklU,t)til 



{hundreds even) 



j=l£=l 
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Eq. (hundredseven) is the discrete analog of (sixteen) and shows precisely how the second 
and other derivatives of the path in that expression are to be interpreted in a time-slicing 
representation. 

We next assume decoherence and carry out the integration over the ^'s neglecting the 
constraints of the coarse graining as discussed in Section III. Note that A in (hundredseven) 
depends on £q om y through the first term and that £jy = because of the 5-function 
in (hundredtwo) . The result for the diagonal elements of the decoherence functional is 
explicitly 



JV-s-00 

N-1N-1 



fc=0 



N-2 N 



*«) = jsL / ( £> ) ( ct ) n «* 






i 



2 



x exp 



\ E E e ' ^(^m)4 nV (j^) e/(X m )]u; (^0,M (^-^)) • 
i=i e=i V V < // 

(hundredeight) 

In this expression det fc/ and /cj nv are the determinant and inverse respectively of the 

(N — 1) x (N — 1) matrix kj(j, £), j, £ — 1, • • • , N — 1. The quantity e^ is the discrete 

version of the equation of motion (nineteen), namely 

— 2±± ^ ^— - J - M,- - J^ e fcij(j, £)X £ (hundrednine) 

and w(X,P) is the Wigner function defined in (tenb). Eq. (hundredeight) contains the 
precise measure for the path integral (eighteen) and the precise meaning of the integral 
over the Wigner function in it. 

(c) Functional Fourier Transforms 

In Section V we utilized a functional Fourier transform of the decoherence functional 
to define a distribution functional for the total force. Here, we offer a more explicit 
definition of what such transforms mean. We consider the case of one-dimensional paths 
for simpilicity. Consider a functional D[£(r)]. On paths that are piecewise linear between 
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time slices to, ti, • • -tjsf = T this defines a function, D(£q, £i, ■ ■ -£n) of the values that £(r) 
assumes on these slices. This function may be Fourier transformed in the following way 

G(R 1 ,---R n ;£ )=Ar J I f[ dtk) ex p(-^ Ylk=l e ^ kRk ) D (£o,£i,"-£aO 

(hundredten) 
where e = t/N is the spacing between the time slices and M is a normalizing factor. The 
inverse of this is 






. TV 

exp ( fc S e ^ fcjRfc ] G ( R i> - - - r n; £o) ■ 

k=\ 



D\i{r)\ = f 5Rexp % - f dt£(t)R(t) 



G[R(t),£q] . (hundredtwelve) 



(hundredeleven) 
It is this latter expression, in the limit where the number of slices, N, goes to infinity that 
gives an explicit meaning to the functional Fourier transform that we wrote in (teng) as 

'0 

In particular, (hundredeleven) defines the measure, 5R. Note that because of the way 
that we have differenced the exponent in (hundredtwelve), £o remains untransformed and 
occurs on both sides of (hundredeleven). 

The normalizing factor J\f is physically arbitrary and will cancel from physical ex- 
pressions such as the definition of the equation of motion in (tenj) and (tenk). However, 
mathematically, it must be chosen carefully in order that expressions like (hundredten) 
exist in the limit iV — > oo. For example, when the decoherence functional is given as in 
(fifteen) and (sixteen) for linear problems the normalizing factor would depend on Kj and 
certainly on e. We assume that such a normalizing factor always exists for interesting 
cases. 

(d) An Operator Expression for the Influence Functional 

The influence functional exp(iW[x'(r), x(r)]/h) was defined by the path integral (three) 
but can be represented as the operator expression (forty three) when the initial density 
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matrix factors as in (six) . We now spell out the details of this connection. By inserting 
complete sets of states of the Hilbert space Hq that are eigenfunctions of the coordinates 
Q , the right hand side of (fortythree) can be written 



Sp 



U T ,o [At)] Pb Ul [x(r)]} = J dQ' f J dQ f J dQ r J dQ 5 [q'j-Qj 



Qf \Ut,Q [ x \ t )] \Qo) \Qo\Uj^q[x(t)]\Q f\ . (hundredthirteen) 

But {Q f\Uxfl[x(T)]\Qo) is just the propagator in the Hilbert space Hq corresponding to 
unitary evolution by the time-dependent Hamiltonian Hn(x(t)) [cf. (fortyone)] over the 
time interval [0, T]. This propagator has an elementary path integral representation whose 
derivation we have reviewed in subsection (b) of this Appendix [cf. (hundred)]. It is 

(Qf \Uto\x(t)]\Qo ) = / 5Qexp(z5'n[x(r), Q{t)]/Ti) (hundredf our teen) 

X ' ' J[Qo,Qf] V } 

where Sq[x{t), Q(t)] is the action (thirtyeight) that corresponds to the Hamiltonian Hq{x{t)). 
The path integral is over paths that start at time t = at Qq and proceed to Qf at time 
T and are other wise unrestricted. Inserting (hundredfourteen) into the right hand side of 
(hundredthirteen), noting [eq (forty two)] that (Qq\pb\Qo) — Pb{Q'qiQq)i an d using (thir- 
tyeight) for Sq[x(t),Q(t)], we recover the expression (three) that defines the influence 
functional exp(iW[x'(r), x(r)]/h). Thus the identity (fortythree) is verified. 
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Figure Caption 

Figure 1. The decoherence of histories coarse-grained by intervals of a distinguished set 
of configuration space coordinates. The decoherence functional for such sets of histories is 
defined by the double path integral of (five) over paths x'(t) and x(t) that are restricted 
by the coarse graining. These path integrals may be thought of [cf. (hundredtwo)] as the 
limits of multiple integrals over the values of x and x on a series of discrete time slices of 
the interval [0, T]. A typical slice at a time when the range of integration is constrained by 
the coarse graining is illustrated. Of course, only one of the distinguished coordinates x a 
and its corresponding x' a can be shown and we have assumed for illustrative purposes that 
the regions defining the coarse-graining correspond to a set of intervals A a , a — 1, 2, 3, • • • 
of this coordinate. On each slice where there is a restriction from the coarse graining the, 
integration over x' and x will be restricted to a single box. For the "off-diagonal" elements 
of the decoherence functional corresponding to distinct histories, that box will be off the 
diagonal (e.g. B) for some slice. For the diagonal elements, corresponding to the same 
histories, the box will be on the diagonal (e.g. A) for all slices. 

If the imaginary part of the influence phase W[x'(t), x(t)] grows as a functional of the 
difference £(t) = x'{t) — x(t), as it does in the oscillator models [cf (ten)], then integrand 
of the decoherence functional will be negligible except when x'(t) is close to x(t) a regime 
illustrated by the shaded band about the diagonal in the figure. When the characteristic 
sizes of the intervals A a are large compared to the width of the band in which the integrand 
is non-zero the off-diagonal elements of the decoherence functional will be negligible because 
integrals over those slices where the histories are distinct is negligible (e.g. over box B). 
That is decoherence of the coarse-grained set of histories. Further concerns the evaluation 
of the diagonal elements of the decoherence functional that give the probabilities of the 
individual histories in decoherent set can be simplified. If the integrations over x' and x 
are transformed to integrations over £ = x' — x and X = (x + x)/2 the restrictions on the 
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range of the ^-integration to one diagonal box may be neglected with negligible error to 
the probability. 
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